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PREFACE 


The theory of tiorbulence has made so much progress during these last 
years that it is of interest to state exactly the obtained results, the 
hypotheses on which these results are based, and the directions in which 
new research is being conducted. 

Messrs. Bass and Agostini have undertaken this work and have summar- 
ized our actual knowledge of turbulence in a series of conferences which 
took place at the Sorbonne, within the Paris Institute of Mechanics. The 
reader will find in the following pages the text of these conferences , 
perfected and revised by Mr. Bass. 

In view of the magnitude of the subject and its simultaneously phy- 
sical and theoretical aspects it had seemed advisable to entrust this 
work to a team formed by a mathematician and a physicist. Mr. Bass had 
taken the responsibility for the theoretical part, Mr. Agostini for the 
physical part. 

Initially, the report was intended to contain three theoretical and 
two physical chapters followed by a chapter on the technique of the 
measurements and on the appratus used: anemometers and statistical- 

measurement apparatus . 

The unexpected death of Mr. Agostini in August 19^9 unfortunately 
made modifications of the original project necessary. This premature 
death deprived us of a highly valuable physicist, at the peak of intel- 
lectual maturity, whose current work on these problems showed particu- 
larly remarkable promise . 

Mr. Agostini had only just begun drawing up the two last chapters; 

Mr. Bass had to take up the editorial work and to complete it according 
to Mr. Agostini’s notes. The chapter on the technique of the measurements 
has been omitted and will form the object of a later publication. 

In order to make up for the gap in the experimental part, the text 
was supplemented by some curves furnished by Mr. Favre which will allow 
utilization for numerical calculations and will enable the reader to 
judge the agreement between theory and tests . 


A. 

Professor 


Fortier, 
at the Sorbonne 
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THE THEORIES OF TURBULENCE* 
By L. Agostini and J. Bass 


INTRODUCTION 


The theory of turbulence reached its full growth at the end of the 
19 th century as a result of the work by Boussinesq (I 877 ) and Reynolds 
( 1895 ). It then underwent a long period of stagnation which ended under 
the impulse given to it by the development of wind tunnels caused by the 
needs of aviation. Numerous researchers, mathematicians, aerodynamicists, 
and meteorologists attempted to put Reynolds' elementary statistical 
theory in a more precise form, to define the fundamental quantities, to 
set up the equations which connect them, and to explain the peculiarities 
of turbulent flows. This second period of the science of turbulence ended 
before the war and had its apotheosis at the 1958 Congress of Applied 
Mechanics . 

During the war, some isolated scientists - von Welzsacker and 
Heisenberg in Germany, Kolmogoroff in Russia, Onsager in the U.S.A. - 
started a program of research which forms the third period. By a system 
of assumptions which make it possible to approach the structure of tur- 
bulence in well-defined limiting conditions quantitatively, they obtained 
a certain number of laws on the correlations and the spectrum. These 
res\ilts, once they became known, caused a spate of new researches, the 
most outstanding of which are those by the team Batchelor-Townsend at 
Cambridge . 

The analysis of these works became the subject of a series of lec- 
tures at the Sorbonne in February -March 19 ^ 9 , which subsequently were 
edited and completed. The mathematical theory of tirrbulence had already 
been published in 19^6 (ref. 5) hut practically Ignored all publications 
later than 19^0. Since the late reports have improved the mathematical 
language of tiirbulence, it was deemed advisable to start with a detailed 
acco\mt of the mathematical methods applicable to turbulence. Inspired 
at first by the work of the French school, above all for the basic prin- 
ciples, then the work of foreigners, above all for the theory of the 
spectrum. 


"Les 

Techniques 


Theories de la Turbulence." Publications Sclentif iques et 
du Minlstere de L'Air, No. 257> 1950. 
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The first chapter deals with the precise and elementary definition 
of turbulence (sections 1 and 2) and describes the tools of mathematical 
statistics on which the ultimate developments are based. Starting from 
paragraph 5^ chapter the reader should be familiar with the defini- 
tions of the calculus of probabilities, the theorem of total probabili- 
ties, ant^ the theorem of compoimd probabilities. This chapter is entirely 
theoretical, and its aim is to review the methods suggested by the theory 
of random functions and which seem likely to be applied to t\arbulence. 

Only the use of Navier's equations has, so far, produced positive results, 
and chapters III, IV, and V are largely devoted to it. However, it should 
be pointed out that there are theories less familiar to hydrodynamiclsts 
which have been proved in other branches of physics (kinetic theory 
of gases, quantum mechanics). The purely random method is described in 
paragraph 5, and its adaptation to moleciilar systems (according to Born 
and Green) in paragraph 7. The statistical method has the advantage of 
furnishing a remarkable demonstration of the general equations of hydro- 
dynamics (paragraph 6) and of providing- an exact classification of the 
statistical parameters of turb\xlence (paragraph k), which is interesting 
to keep in mind when studying the foreign reports, too exclusively devoted 
to spatial correlations. In any case, the reader who wants to read 
chapters IV and V can pass up most of chapter I, except perhaps para- 
graphs 1 and 2, without major trouble. 

Chapter II deals with the kinematics of statistical medi\jms and, 
particularly. Isotropic mediums. It seemed practical to include at the 
same time the velocity correlations, the theory of which has been given 
in almost final form by Karman, at the end of the second period, and of 
the spectrum, the theory of which, due to Taylor's initiative, has only 
been achieved very recently. Only paragraphs 10 and 11 refer to statis- 
tical functions, and their detailed knowledge is not indispensable for 
reading the rest of the chapter. The results and the formulas of 
chapter II are constantly applied in the subsequent chapters, but it is 
not necessary to know the proofs which are, in most cases, a simple matter 
of calculation. The most important of these formulas are, moreover, com- 
piled in a special section following chapter V. 

Chapter III is a mathematical study of the application of Navier's 
equations to turbilLent motion. Paragraphs 17, I8, 19, and 22 are f\mda- 
mental. Their main purpose is to recall Karman ’s results of 1938 with 
some improvements and some supplements of more recent date . The para- 
graphs 20 and 21 review a certain number of physically reasonable solu- 
tions of the Karman-Howarth fundamental equation. Some of these assume 
particular importance in chapter V but, first, it seems advisable to give 
an impartial view of the whole and to proceed progressively into the^ 
domain of the concrete. The paragraphs 25 and 2k deal with an equation 
by Heisenberg which involves time correlations and from which probably 
not all possible results have been extracted. It is not indispensable 
to have knowledge of this in order to continue. Paragraph 25 contains a 
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mathematical accoimt of Heisenberg *s numerical theory of the spectrum^ 
which is taken up again in chapter V from a more physical point of view 
and whose examination^ contrary to paragraphs 25 and 2k, proved useful 
before attacking chapter V. 

Chapters IV and V deal with new physical theories involving simi- 
larity hypotheses and producing numerical laws. Chapter IV reviews the 
works of Kolmogoroff and Weizsacker^ chapter V those of Heisenberg^ 
Batchelor^ and Townsend on the decay of turbulence created by grids. 

Finally, in an appendix, the theoretical discussions of chapters III, 
IV, and V are illustrated by some correlation curves and spectrum curves 
measured directly in the wind tunnel by A. Favre, in the laboratory of 
the mechanics of the atmosphere at Marseille, or derived from experi- 
mental curves by elementary transformations. 

An exhaustive study of modern theories of turbiilence calls for some 
knowledge of the calculation of tensors, probabilities, and statistical 
analysis, besides the classical conceptions of differential and integral 
calculus . 

As regards the tensors, knowledge of the definitions and fundamental 
operations with rectangular cartesian coordinates is sufficient. There 
are a number of articles on this subject, but they generally lean toward 
the tensor analysis with curvilinear coordinates for which there is no 
need, (it should be noted that the tensor analysis plays, in contrast, 
an important part in the theory of the boundary layer around an airfoil. ) 
Incidentally, there is available a little book recently published, by 
Lichnerowicz, entitled: Elements of tensor calculus (collection 

Armand Colin). On mathematical statistics, the book by Darmois, pub- 
lished by Doin (1928), can be consulted. For the elementary theory of 
random functions, consult the first part of Bass* report (ref. 5) 
the appendix to d*Angot*s "complements of mathematics" (editions of the 
Revue d*Optique), edited by Blanc -Lapierre . More detailed information 
on statistical functions can be found in Levy*s book: "Stochastic 

Processes and Brownian Motion" (19^8). 

The present report deals only with general theories which are valid 
whatever the physical or geometric causes of turbulence may be, as is 
shown in chapter IV. Only in chapter V the assumptions are limited and 
the study involves a problem of decay of turbulence that is compatible 
with the turbulence in wind tunnels. These theories are probably appli- 
cable to mediums of extremely diverse scales, from the microturbulence 
to the terrestrial atmosphere (Dedebant and Wehrle), to stellar atmos- 
pheres and to interstellar matter (Weizsacker ) . However, among the 
hypotheses there is always that of the incompressibility, which precludes 
the application to sonic or supersonic flows. This is not an indispensi- 
ble hypothesis, but it simplifies the calculations considerably, and the 
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consequences are easy to check. So the incompressible turbiilent motions 
are practically the only ones studied up to now. 

Among the other hypotheses worthy of discussion, those referring 
to processes of energy transport (paragraph 2 h and chapter V) assume that 
all turbvLlent energy comes from the motion of the whole, dissociated an d 
broken up by the obstacles of periodic structure. But nature furnishes 
examples of different turbulences such as the so-called thermal turbu- 
lence, for example, where the source of energy is not the fain of a wind 
tunnel but the solar radiation suitably transformed into kinetic energy. 
Kolmogorof f ' s theory of local isotropy applies probably to thermal turbu- 
lence, but the forms of energy transport in the spectrum must be differ- 
ent from those encountered in wind tunnels. The same applies to the 
astronomic turb\ilent medl\ams alluded to previously. So the foregoing 
remarks limit the scope of the recent theories, in a certain measure. 

Omitted entirely was the problem of turbulent boundary layer which, 
experimentally, depends on the same technique, but has been approached by 
different mathematical methods. Furthermore, it is a complicated problem 
where the turbulence is neither homogeneous nor isotropic . 

In what measure are the results, suitably demarcated, defined? The 
statistical character of the tiirbulent velocity seems a clear and well- 
established notion and consequently the mechanics of turbulence will be 
a statistical mechanics. As far as the kinematics (chapter II) are con- 
cerned, we are therefore on solid ground. The dynamics of turbulence 
(chapter III) itself is likewise well established, by means of the 
hypothesis of the validity of the Navier equations. This hypothesis, 
generally adopted because it is convenient and, one might say, necessary, 
has however at times raised considerable doubts. However, while experi- 
mental verifications do not contradict it, there is yet no occasion to 
reject it. 

But, what should be expected from experimental verifications? First 
of all, it is fo\md that the theories are still rather imperfect. In 
fact, the theories are usually limited, acceptable in the limiting con- 
ditions which are difficult to attain actually (very high Reynolds num- 
bers, for example). No rigorous verification should be expected since 
the true conditions are too far removed from the theoretical conditions. 

On the other hand, the accuracy of measurement is low. The original 
reason for it lies in the very nature of the turbulent phenomenon, and 
its irregvilar and badly defined character. Fiarthermore , the anemometers 
are coarse Instruments, their operation not s\off Iciently kn own in the 
presence of turb\jlence, and their interference with the fluid sometimes 
a little mysterious. All this helps to lower the experimental precision. 
To measure the "length of dissipation," for example, two ways are open; 
one is to measure the correlations of the velocity at two infinitely close 
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points which has no practical sense since two anemometers cannot he 
brought together indefinitely; the other is to use one anemometer for 
measuring the mean square value of the derivative of velocity. But the 
operation of taking a derivative of a function as complicated as the 
velocity is by itself inaccurate and introduces serious scattering of 
the points of the derivative ciorve. So in both cases, the accuracy of 
the results is extremely limited. 

In consequence, the experimental verifications can be applied only 
to the orders of magnitude. To illustrate; to check whether a parameter, 
according to theory, is constant, one constructs the representative curve 
of the parameter and, if this curve has a sufficiently extended maximum, 
one estimates that the experiment closely confirms the theory, conceding 
that, as the experimental conditions more nearly approach those stipu- 
lated by theory, the maximum flattens out more and more, and one does 
not appear too severe in the examination of the scatter of the test 
points . 

On this assumption, the verifications of the experimental laws of 
turbulence, due in particular to Townsend, are encoiaraging. Therefore, 
it is well to retain the hypotheses of chapters III, IV, and V, although 
some of them obviously have their limitations, and we o\ight not hope for 
results greater than actually can be given. Take an example drawn from 
the theory of the spectrum, for instance. At the beginning of decay of 
turbulence in a tunnel and at high Reynolds numbers, a certain spectral 
function F(k) is of the form Ck^ for the small values of wave num- 
ber k. , It then passes through a maximum, then becomes proportional 
to k”^' Sind finally approaches zero as k — > », maybe as k"”^. But 
the regions of the axis of k in which the fragments of the laws to be 
emmierated remain acceptable are badly defined and connected by zones 
of which the structiire is not known. So the future task of the theorists 
will be to combine these partial res\ilts into a single acceptable law at 
least to the extent that it does not become fundamentally incompatible 
with the nature of turbulence, since k = 0 up to k = ». Only then 
wi3J. there be a true theory of turb\ilence. 
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CHAPTER I 

GENERAL ASPECTS OF TURBULENCE 
THE STATISTICAL METHOD 


1. Definition of turbulence: 

Our purpose is to discuss the operation which consists of measuring, 
at a given point, the velocity of a flow which, to avoid every difficulty, 
is assumed to be steady and uniform at the \isual macroscopic scale. The 
measioring Instrument is an "anemometer" with approximately determined 
dimensions and time constaint, and of sufficiently ideal nature so as not 
to distijrb the flow by its presence. If this anemometer is large enough, 
it measures the "velocity of the main flow" of the fluid. If it is very 
small, it can be imagined that it operates in discontinuous manner, never 
undergoing the influence of more than one molecule at a time . An anemom- 
eter sensitive to the individual action of molecules is, of co\irse, 
unattainable, but the idea of such an instilment is convenient for repre- 
senting the extreme limit of fineness of kinematic measurements in a fluid. 
Between these two extremes, the indications of the anemometer depend upon 
the structure of the fluid. It may happen, by exception, that, when its 
dimensions are progressively reduced, the velocity which it indicates 
remains unchanged up to the moment where the individual influence of the 
molecules starts to mak e ItselP felt and where the indications lose all 
statistical significance. The flow is then said to be laminar. 

But, in general, the matters are otherwise. We start with a first 
anemometer which, through its dimensions, fixes a certain scale of meas- 
urement. This anemometer measures the mean speed of the molecules in a 
certain volume V. Then it is replaced successively by smaller anemometers 
in such a way that volume V decreases progressively. It happens that, 
from a certain value V]_ of V, the numerical indication supplied by the 
anemometer changes. If V is decreased continuously, the new indication 
remains stable up to a certain value V2^ then changes again and so on. 

The intervals (V^, V2), (V2, V5), . . . characterize the various scales 

of turbulence, and the motion of the fluid is said to be turbulent. More 
exactly, they are the conditions necessary for a fluid to be turbulent, 
and which must be defined and perfected to make them sufficient and 
practical. 

The last value Vn of the series V^, V2, ... is that from which 

onward, the notion of average loses its significance, the number of mole- 
cules contained in the volume not being large enough any longer for 

statistical purposes. The series V2^ • • • can be discrete or 
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continuous. If it is discrete^ it still does not imply that the critical 
values V 2 , • • • are mathematically defined. In the vicinity of 

a rapid variation of the anemometer indications occurs, which subsequently 
become stabilized in the region ^ 2 )^ and so forth. If this stabili- 

zation is not very clear, the turbulent scales are said to succeed one 
another in continuous fashion. 

The previous discussion ends with the notion of the turbulent fluid. 
But the definition of the laminar motion given above is a little too 
restrictive and the distinction between laminar and turbulent still not 
precise enough, as proved by the following example: 

Consider the motion of air produced by stationary waves in a sound 
tube. The motion of the whole reduces, at rest, to large scale. But, 
at each point there exists a speed other than zero, a periodic time func- 
tion which, at a given instant, varies periodically from one point to 
another. This motion, lying between the system at rest and the molecular 
agitation, has not a turbulent character. 

Turbulence, as shown, implies first the notion of scale. But it 
should be added that, at a given scale, each component of the velocity 
at a point is a function of time presenting a character of periodicity 
without fundamental period. This is not a periodic function but a sum 
of harmonics, the frequencies of which are not multiples of an identical 
fundamental frequency. This irregularity of the turbulent agitation is 
essential and distinguishes it from sound agitation, or preturbulent 
vortical motions, like the cellular vortices of Benard. The mathematical 
symbol for the turbulent velocity is not the ordinary Fourier series, but 
Fourier *s integral. It will be discussed later. 

This concept of irregular agitation at a point as function of the 
time is not itself sufficient. It makes it possible to differentiate the 
turbulent agitation from the periodic sound agitation (musical sound), 
but not from the noise, which is an agitation without definite period. 

What distinguishes the noise from turbulence is the fact that it is prop- 
agated by waves that exist on surfaces of equal phase, and consequently 
have a regular spatial distribution, notwithstanding the irregularity in 

the time of the local velocity.^ 


^A descriptive and purely kinematic distinction is involved here. 

Its cause (compressibility) is not discussed. 

The difference between turbulence and sound agitation should become 
plain from the following example: In a turbulent wind tunnel, we select 

at a point the longitudinal component of the velocity with a hot wire and 
send the electric current of the hot Wire to a loud speaker on the outside 
of the tunnel. The atmosphere becomes the source of an irregular agi- 
tation, which propagates by waves and is not turbulence, although, at 
a point, the internal turbulent motion and the external sound motion 
have some important kinematic elements in common. 
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For a given scale, each component of the turbulent velocity is an 
irregularly periodic fimction of both space and time. 

It seems that turbulence is well defined by these kinematic condi- 
tions and hence is distinguished from all other more organized fluid 
motions . 

2. Average values - statistics: 

Figure 1 represents the record of a component u of the turbulent 
velocity as function of the time. (Record of the velocity of turbulent 
agitation in a 20-cm by 50 -cm tunnel: airspeed, 20 m/sec; intensity of 

turbulence, 5-10”^; time of recording is 0.05 second.) 

The most natural method of measuring the mean velocity on the graph 
consists in forming the ordinary integral 


u 


m 



u(t) dt 


( 2 - 1 ) 


extended over the total duration T of recording. This method is, in 
general, not very satisfactory because the operation lacks precision 
when the curve u(t) is complicated. 

A more precise method consists in dividing the graph by parallels 
to the axis of t, suitably close together in the ordinates U 2 ^, U 2 , 

Uj, . . . , in measuring the number nj^ of points where the line of the 
ordinate + Ui+]_) meets the curve u(t) and then in computing the 

quantity 

u = ^ ui ^ (2-2) 

1 ^ 


n = XZ is the total number of points met by all parallels, 

i 

For this calculation, a profitable first stage consists in first 
constructing the graph giving the corresponding statistical frequency 

= Zi for each velocity Uj^. Crossing the limit obviously makes it 

possilhLe to plot a curve of frequency f = f(u) (fig. 2) such that the 
proportion of the values of the velocity comprised between u and u + du 
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is equal to f(u)du, the integral 


0 


f(u)du, which replaces 


being equal to unity. The final expression of u is then 



u 


p 00 

/ uf(u)du 


(2-3) 


It replaces (2-2) and shovild be compared with (2-l). 

The practical operations enabling the replacement of by u 

correspond to well-known mathematical operations. The mode of com- 
puting u is that of a Lebesgue integral, and Um is an integral of the 
classical type of Riemann. If u^ is computable, both methods yield 
the same result. But it may happen that Riemann 's integral does not 
exist because the function u(t) is too complicated mathematically. 
However, in general, Lebesgue 's integral exists (if u(t) is measurable, 
and natirrally boimded). This mathematical case corresponds to the 
practical case where the curve u(t) is too complicated for an accurate 
continuation of the integration. The function f(u), continuous and 
differentiable in the cirrrent cases, is the medium which, determined once 
for all, replaces the calculation of Lebesgue 's integral by that of an 
ordinary integral by means of the plotting of curve f(u). 

From the function f (u), other averages can be computed. For 
example, the amount of differences of the speed with respect to_its 
mean value can be figured by computing the mean value of (u - u)^. 

Rather than defining this average by Riemann ' s integral 


T Jq 

[u(t) - uj^dt 

(2-4) 

more precise to 

use the formula 


(u - u)^ = 

r“ 2 

/ (u - u) f(u)du 

' —oo 

(2-5) 


which, once the curve f(u) is plotted, calls only for operations of a 
simple character. 
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Thiis, it is seen how much the construction of the curve f (u) 
simplifies the numerical calculations of turbulence. Various experi- 
mental techniques make for direct attainment of this curve without passing 
through the numerical analysis of a velocity record. For measuring the 

averages such as (u - u) , it is often advisable also to use specialized 
equipment without first plotting the curve f(u). 

5. Random variables and the laws of probability: 

The theoretical significance of the function f(u) is analyzed. 

It groups the statistical data contained in the initial curve u(t), with 
this exception that the chronological order in which the velocities actu- 
ally follow one another does no longer appear. This limitation is quite 
natural, though, and it will be seen later that this order reappears, in 
a certain measure, by the introduction of space and time correlations. 

Obviously, only statistical data can supply stable information on 
turbulence. When the same record of the velocity is begun again several 
times while taking every reasonable precaution so that the conditions 
are identical, it obviously results in curves u(t) which absolutely 
are not superposable . The function u(t) has not, therefore, the char- 
acter of permanence that is suitable for representing the laws of a 
physical phenomenon. But this character is relevant to the function f(u) 

or to averages such as u, (u - u) whose values are characteristic 
nvimbers of the investigated flow, and which are derived by simple mathe- 
matical operations from the f\mctlon f(u). Hence, we direct our atten- 
tion to this function which can be regarded as representing the first 
law of t\irb\ilence . 

To say that the velocity is characterized by a curve of statistical 
frequencies is to say, by comparing the frequencies with probabilities 
nnH f(u) with a density of probability, that this velocity is a statis- 
tical quantity; f(u)du is the probability that the chosen component of 
the velocity is contained between u and u + du. 

A priori , such a law of probability could be dependent on the time . 
That would correspond to a turbulent flow for which the laws varied with 
respect to time. There is no contradiction to the initial assumptions 
of permanence here. It is a question of scale. In order for the experi- 
mental operation by which f(u) is defined to have any meaning, it is 
necessary that two conditions be realized sim^lltaneously: 

(1) The number of oscillations in the time interval T involved 
must be great enough to furnish satisfactory statistics. 

(2) The laws of turbulence in this time interval T 
ally permanent. 


must be practic- 
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If it is not so, it might he difficult to reconcile the statistical 
theory with the experiment. Fortunately, those "ergodic" conditions are 
practically always realized in the usual cases and are therefore taken 
for granted in the following: 

The velocity has three components uj_, U2, and u^ which are treated 
as three random variables, components of a random vector. 

It should be noted here that the velocity is perhaps not sufficient 
for characterizing the turbulence. It might appear useful to Introduce 
other quantities, such as the pressiore, which should be treated as a 
statistical quantity. But, owing to the equations of motion, this then 
will be a function defined by the velocity and its derivatives. For the 
present, it is assumed that the txirbulent motion is sufficiently well 
defined by its velocity so that the problem narrows down to the laws of 
probability applied to the velocity. 

The simplest of these laws, that which immediately generalizes the 
experimental function f(u), is the law of probability of the system of 
three components u^^, U2, u^ of the velocity. This law may vary as 

function of the time t (problem of spontaneous decay of txirbulence) and 
of the space (variation of turbulence in terms of the distance from the 
walls). It is therefore a function of t and the ordinates xj^, X2, xj 
of the point of measurement. Its density is denoted by 

f^ui, U2, ny, x^, X2, Xj, t) 

or, abbreviated^ t). 

The quantity f du]_ du2 du^ or, abbreviated, f du, represents the 
probability that, at the point xj_, X2^ (or x) and at the instant t, 

the three velocity components are comprised between and + duq, 

U2 and U2 + ^^2^ ^5 du^. 

But a single law of probability defined in terms of four parameters 
Xj_, X2; x^, and t is not adequate for characterizing turbulence. It 
is necessary to introduce the more profound concept of random function 
and to consider the turbulent velocity as a random function of space and 
time. This is the random velocity field. 

This point is now to be defined. An isolated random quantity is 
defined by its law of probability. But, to define a system of coexistant 
random quantities requires more than just their laws of individual prob- 
abilities. The stochastic dependencies or correlations between these 
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quantities must also be known. Taking a family of probability laws does 
not give the right to speak of the system of corresponding random quan- 
tities without completing the data by those of the correlations. Sup- 
posing that the fajnily in question depends on a continuous parameter t. 

We know then an isolated random quantity U(t) for each value of t. 

This immediately suggests grouping the U(t) corresponding to the var- 
ious values of t in a well-defined system of statistical quantities. 

This calls for the introduction of the correlations between the U(tq), 
U(t 2 ), . . . corresponding to an arbitrary system E of the values tq, 
t 2 , ... of the parameter. To proceed thus, means to define a random 
function. Natxirally, this also holds for laws of probabilities dependent 
on several parameters. Thus, when a turbulent medium is represented as 
a velocity field, the system of velocities at each point and at each 
instant precisely constitutes a system of coexistent statistical quantities, 
of which the physical interactions characterizing the structure of txorbu- 
lence have the correlations for mathematical description. 

Among the systems E, the simplest are the denumerable systems and 
even the finite systems and, among the latter, the simplest one which is 
not trivial is that of two elements, that is, of two points of space and 
time . 


The concept of random fimction thus sviggests the comparison of the 
velocity vectors at two different points of their field of definition, 
that is, for two different positions x and x' (x represents the 
point of the coordinates Xq, X 2 > x^) and for two instants t and t'. 

It concerns a statistical comparison which makes it possible to define 
the law of probability of the two systems of the velocity components at 
points X, x' and instants t and t'. This law has a very clear 
physical meauiing and is easy to define experimentally or, at the least, 
to construct the surface of (statistical) frequencies corresponding to 
one velocity component at point x and a second component at point x', 
the measiirements being spaced at a chosen time interval r . 

In practice, the question is frequently handled from a less general 
point of view. One is not concerned with the laws of probability them- 
selves as such but only with their most simple moments, those of the 
second order which are associated, as will be shown, with certain 
"physical" aspects of turbulence and, possibly, with certain moments of 
the third order that play a part in modem theories. The moments of the 
second order constitute the correlation tensor of the law of probability 
of the velocity field at two points and at two instants. They are meas- 
ured direct, without resorting to frequency curves or frequency surfaces 
or to velocity recordings. A detailed study follows later. 
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4 . The concept of random point - velocity field - turbulent diffusion: 

The density of probability f(u; x, t) of the velocity field contains 
the velocity of the whole fluids but not its density. With p(x3_, X2, x^; t 
or^ abbreviated^ p(x, t) denoting the quotient of the density at point x 
by the fluid mass^ p is a normalized function^ as a density of probabil- 
ity^ which means that 


where dx represents the element of volume dxj^ ^2 ^^5 integral 

is extended to the volume V occupied by the fluid. 

In all modern studies on turbulence, the fluid is naturally assumed 
incompressible, so that p is a constant, equal to l/v in the volume V 
and zero at the outside, p could be simply replaced by a constant; but 
the more general conclusions to be arrived at ultimately are more complete 
if this simplification is not made. On the other hand, it is interesting 
to foresee, at a certain stage of the theory, the day when it will be 
possible to study the turbulent motions in conditions where the compressi- 
bility is no longer negligible. For these reasons, p is treated here 
as a function of x (and even of t, if necessary). 

The product R(x, u; t) = pf is now formed. It obviously is 
normalized with respect to the system of the six variables x, u 


R presents thus the characters of a density of probability with respect 
to these six variables. The quantities Xj_, X2^ x^ are regarded as the 
coordinates of a moving point, U3_, U2, as the velocity components of 

this point. These are six random quantities of which the law of probability 
at instant t is known. Thus a random point can be associated with the 
turbulent fluid in correspondence with a given scale. 

This point is now to be discussed as was the velocity field in the 
preceding paragraph. 

The position and the velocity of this point are random functions of 
the time. The theory of random functions suggests the study of the law 
of probability of the system of positions and velocities of this point 
for an arbitral^ combination of instants. The initial analysis was on 


p dx = 1 


(4-1) 


V 



(4-2) 
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the law of probability at one instant. The generalization from one to 
two instants seems to us adequate for forming a physical theory of t\rrbu- 
lence and, in particular, for considering the statistical organization 
in time of the velocity, and what may be called the "interactions between 
txirbulent particles." 


The values of the positions and of the velocity must therefore be 
associated to two instants t and t'. This association is a stochastic 
relationship defined by the density of probability of the following system 
of 12 statistical variables: 


x^, X2, 

"5 

position 

Uj_, u^. 

^5 

velocity 

Xj_', X2' 

' X ' 

position 

Ui', U2' 

' 11 ' 

velocity 


> at instant t 


i 


at instant t * 


or^ abbreviated^ 


G(x, X* , u' ; t, t * ) 


and it is assumed that it defines the turbulent motion. By what stages 
caji it be measured ? 

According to the theorem of compound probabilities^ G dx dx* du du* 
is the product of two factors: 

(1) The probability that at the instants t and t* the statistical 
image point of the fluid might have positions contained within the inter- 
vals (x, X + dx) and (x*, x* + dx*). 

(2) The probability (conditional) that^ these positions being fixed, 
the velocities are contained in the intervals (u, u + du) and (u*, u* + du'). 

This last probability is designated by H(u, u*; x, x*, t, t*)du du*. 
According to the theorem of compound probabilities, the probability (l) 
is the product of the probability p(x, t)dx that the random point is 
found, at instant t, in the interval x, x + dx, through the (conditional) 
probability that, the position at instant t being chosen, its position 
at instant t* is found in the interval x*, x* + dx*. 


NACA m 1377 


15 


This last probability is designated by p(x'j t, t')dx'. We can 

write 


G(x, x', u, u', t, t') = H(u, u'; X, x', t, t')p(x'; x, t, t')p(x; t) 

(^- 3 ) 

What is the significance of the three factors of which G is the 
product ? 

We already know p, which represents except for a numerical factor 
the density of the fluid. 

The function H is the law of a random velocity field at two points 
and at two instants. It is natural to identify it with the law of the 
field already discussed at the end of paragraph 5* It is seen that it 
gives no complete picture of turbulence. It does not explain the func- 
tion p. 


The function p represents the turbulent diffusion at a chosen 
scale. It is the relative density at point x' and at instant t* of 
fluid elements which have passed neighboring point x at instant t. 

It should be pointed out that the diffused portion of the fluid is essen- 
tially compressible since the density p(x^; x, t, t^) decreases in 
proportion as the point x* is removed from the initial point x where 
it is maximum. The spread can be materialized and p can be measxired 
by introducing with the necessary precautions at point x a dye that 
spreads in the fluid. But it is a rather ticklish matter to separate 
the effects of the various turbulent scales^ especially of the molecular 
diffusion. It is accomplished by adapting the particles of the "dye" to 
the chosen scale. This way Kampe de Feriet rendered the turbulent dif- 
fusion of the tunnel flow visible in his experiments at the Institute of 
Fluid Mechanics ; at Lille ^ by injecting soap bubbles at a point in place 
of dyed particles. These soap bubbles^ because of their size^ were 
sensitive to the turbulent fluctuations and were used successfully for 
measuring a density of turbulent diffusion. 

5. Equations of development of the laws of probability: 

The analysis of the turbulent velocity field made it possible to 
represent a turbulent fluid by two random functions X(t) and U(t) 
playing the part of the position and of the velocity for a random material 
point. It was shown how this concept of random point gives a very com- 
plete picture of the turbulence. But this picture is qualitative and 
must be made more quantitative. 
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At the heginning, no distinction is made between the velocity 
vector X(t) and the position vector U(t). The six components of these 
two vectors are considered as those of a unique vector X(t) in a six- 
dimensional space, by putting Uj_ = xj^, U 2 = x^, u^ = Xg, and designating 
the density of probability of this vector by R(x; t), this notation being 
the abbreviation of X 2 ^ x^, x]^, x^, xg; t), that is, of 

R(x^, X2, Xj, v . 2 , u^; t). 

We already had applied (4-5)^ the theorem of compound probabilities, 
to the law of probability G of the position and the velocity at two 
instants t, t': 


G(x, x', u, u'; t, V) = jp(x; t)p(x'; x, t, t'^H(u, u'; x, x', t, t') 

The question involved essentially the separation of the velocity, 
which figures in the factor H, from the position. But there is another 
way of applying this theorem. It consists in separating the two 
instants t and t* by writing 

G(x, x^, u, u*; t, t') = R(x', u*; t*)K(x, u; x*, u*, t, t') 


(5-1) 

hence it results, according to the theonn of the total probabilities, 
that 


R(x, u; t) = jR(x', u'; t')K(x, u; x', u', t, t')dx' du' 

(5-2) 

K is a conditional density of probability, that of the "probability of 
passage" of the state of the random point at instant t' to its state at 
instant t > t ' . 

Abbreviated, we get 


R(x; t) = t')K(x; x', t, t')dx' (5-5) 

To exploit this equation, recoiirse is had to a method patterned after 
the concept of J. Moyal (ref. 53)- This idea was to consider (5-5) as a 
linear integral transformation for passing from the density of proba- 
bility R at instant t' to the same density at instant t. The density 
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of probability of the passage K is the kernel of the transformation. 
With K-tt' indicating the linear operation which has as kernel the 
function K(x; x', t, t'), one may write symbolically: 


R(x; t) = R(x’; t') 


( 5 - 4 ) 


This transformation has special properties rather difficult to 
define. Suffice it to state that it reduces to the identical transforma- 
tion for t = t'. It is assumed that it has an inverse 
without, however, prejudicing the relations between K-^t* and 

The infinitesimal transformation applied to 1 is now examined. 

To this end, R is assumed differentiable with respect to t and Sr/^ 
calculated 


dR(x, t) 
^t 


lim 


R(x; t + t) - R(x; t ) 


T — > 0 


T 


= llm ^'t 

T — ^ 0 ^ 


= lim ^'t - 

T ^-0 


If, as assumed, the function R is differentiable with respect 
to t, the limit of the second member exists, and it is a function of x,t 
independent of t*. It can be computed by giving parameter t* (inde- 
pendent of variable t) any value not exceeding t, such as t* = t, for 
example . Hence 


St 


lim 

0 


^t+rt 

T 



( 5 - 5 ) 


The operator 


L = 


'<U.Tt - 1 


lim 
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defines the infinitesimal transformation of K, sind R confirms the 
fundamental functional equation 


which, theoretically, enables R(x; t) to be computed at instant t 
when R(x; t') is known at an initial instant t'. 

The foregoing calculation cannot be explained in a simple manner 
from formula (5-5) because the limit of the function K(x, x', t, t') 

does not exist when t' > t; this is a symbolical "function of Dirac" 

which expresses the identity =1 in the functional formalism. It 

is preferable to pass, as Moyal did, from the densities of probability 
to characteristic functions. Moyal 's calc\ilation follows: 

With the function K is associated the characteristic function of 
the increment x(t) - X(t'), the value of X(t') once fixed, or by 
definition 



(5-6) 


It (a; x' , t, t ' ) = 



X', t, t')dx (5-7) 


One assimies likewise: 



the characteristic function of X(t) . 


According to the theorem of total probabilities , one has : 



x' , t, t ' )R( x' ; t' )dx 


(5-9) 


This relation replaces (5-5) • 

The two members of (5-9) are differentiated with respect to t 


lim ie 
8t T-^0^ 


iox' o(n^ x', t + T, t') - o(oLj X, t , t ) 


R(x' ; V )dx 


T 
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Whereas the conditions of differentiability under the J' sign 

^permutation of signs lim and of J''^ cannot be verified on (5-3), they 

can be here, in all current cases. Since the limit must not depend on t', 
let t' = t, so that 

rei<^R(x', t) lim X', t + t, t) - 1 , s 

at J ^^0 T 


because il>(a; x, t, t) =1, for t' = t. 

The limit that figures under the J' sign is a certain function 
9 (a; x', t'), whence follows the basic equation 


^qp r iox' 

dt J 


R(x', t)0(a; x', t)dx' 


equivalent to (5-6). From that, two probability densities can be recovered 
by taking the Fourier transforms of the two members. Lastly, dR/dt is 
expressed in form of an integral transformation of R, equivalent to the 
transformation L, which, after introducing a kernel function l(x, x', t), 
gives 


6R(x, t) 
dt 



t)L(x, x', t)dx' 


(5-11) 


The second member of (5-11) is represented in integral form but, in 
many cases, it can be expressed in form of a differential operator of 
finite or infinite order. 


Examples . - Supposing the probability of passage obeys the Laplace- 
Gauss law and, to avoid any confusion of the notations, tQ now denotes 
the previo\is instant; the differences of X^tQ^ and X(t) are denoted 
by Sq = S ^tQ^ and S = S(t), the correlation coefficient between X^to^ 
and X(‘t) ^ It is known that the law of probability 

related to X(t), when the value Xq of x(to) given, has for density 


K 


(x; 




to, t 


/irtS\/l - r^ 


2 (l-r^) 


— - r- 


.^0 

’^0 


(5-12) 
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Using Moyal's method, it is shown that the prohahility density R(x, t) 
which assumes K as kernel of the prohahility of passage satisfies the 
equation of partial derivatives 


St 



(5-15) 


where r’ represents the value of — ^ and S* 

Stp ' ^ 

denotes the derivative — S(t) . 

dt 


The equation (5-15) is, moreover, demonstrated very simply by a 
direct method. By definition of K 


R(x, t) 



(5-14) 


this equation being, in particular, satisfied when 


R(x, t) 



When the two memhers of (5-14) are differentiated with respect to 
t and X, equation (5-15) is verified. Reciprocally, the integral of 
(5- 15 ) which is reduced to R^x, tQ^ for t = tQ, can he put in the 

form (5-14) . 

The construction of the random functions compatible with these laws 
of probability is an easy matter. 

To illustrate: 

Let h(s) be a variable random function of the parameter s, obeying 
the reduced Laplace-Gauss law (h = 0, 1? = l), and so that the increments 
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^ (®l) ~ ^ (®2) (^ 5 ) 

(si, S2j and JS5, 34^ , 


- h ^34^ , corre 3 ponding to two 3 eparate interval 3 
are independent. The random function 


X(t) = 



(t - s)dh(s) 


(5-15) 


is now considered. 

It is easily shown^ that X(t) 


for typical difference S(t) = 


t5/2 

VT~ 


obeys a Laplace -Gauss law having 
and that, if t = t - tQ 


r 





(5-16) 


Without going into details, it is simply recalled that the demonstra- 
tion utilizes as intermediary the characteristic function of X(t), of 
which the logarithm, owing to the properties of independence of the 
dh(s), is expressed by an elementary integral. 

From these formulas, it follows that 


r* 


= 0, 


SI = 5 1 
S 2 t 


and that R(x, t) 


verifies the first-order partial differential equation 


M + 5 1 ^ 0 

ht 2 t Sx 


^See, for instance, Bass, reference or ”The random functions and 
their mechanical interpretation." Revue Scientifique, No. 52^0^ 19^5- 
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X(t) is a differentiable random function. Its stochastic deriva- 
tive is deduced from the expression (5- 15) operations of classical 

r ^ 3 X 

form and written as / dh(s). It is easily proved that — — is the 

J 0 2 t 

related mean of this derivative, if the value x of X is fixed. The 
equation (5-17) bas many solutions which are densities of probability. 
It does not determine the function K. 

Next consider the elementaiy random function 



(5-18) 


stochastic derivative of the function which has been studied as the first 
illustrative example. It is shown that X(t) obeys a Laplace-Gauss law 

having S(t) = /b for typical difference and that, if tg < t: 


r 




^ (5-19) 


Therefore, 

equation 


R(x, t) verifies the second-order partial derivative- 


M i ^ 

St 2 Sx^ 


( 5 - 20 ) 


This is the equation of heat. The solution which reduces to a 
given function Rq(x), for t = t^, is given by the classical formula 


R(x, t) 



which can also be shown by using the Fourier transform of R, that is, 
its chajractistic function. 
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Thus equation (5-20) defines here the form of the probability of 
passage K^x; Xq, t, tQj contrary to what occurred in equation (5-17). 

This example presents an unusual peculiarity. It is easily verified 
that the function 


K 


(x; 


O' 





(5-21) 


satisfies the functional equation 


p 00 

K(x 2 ? Xq, t2, to) = / K(xp; Xi, tp, ti)K^xi; xq, t]_, to)dX]^ (5-22) 

O _ 00 


called the Chapman-Kolmogoroff equation, which characterizes the 
Markoff processes (or more generally, the ”pseudo-markovian" processes), 
a functional generalization of simple Markoff chains. This equation 
can be written in operational notation as 


^tpto ^ 2^1 ^''"l‘'"0 


(5-25) 


It expresses that the operations K form a group. The operator 


^+Tt' ^'t ~ ^ 
T 


is then written simply as 


Kt+Tt - 1 

T 


It is independ- 


ent of t ' . In the general case it is its limit L only when t — > 0 
which must be independent of t, but it itself is not. 


Returning to equations (5-12) and (5-15) it now is assvimed that 
R is not dependent on t; S is then a constant and S' = 0. If 
r' ^0, eq\iation (5-15) is written as 


A(xR) + S2^ = 0 

hx ^x2 


(5-24) 
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The only solution of this differential equation which defines a 
law of probability is the density of the Laplace -Gauss law 


R(x) 



(5-25) 


This example contains, as a special case^ stationary random functions 
for which the function r^t, tQ ) depends solely on the difference 
t - tQ; r' is then a constant, independent of t. 

The final example deals with a vectorial random function of a type 
to be utilized later. Consider simultaneously the random function 



(t - s)dh(s) 


and its derivative 


X'(t) 



(5-26) 


( 5 - 27 ) 


which can play the part of the speed. 

It is easy to form the 'functions R and K for the vector X 
having X and X’ as components. To find the partial differential 
equation verified by the function R(x, x^, t) it is not necessary to 
first form the kernel K. It is simpler to begin with the expressions 
of X and X^, which, passing through the intermediary of the charac- 
teristic function of X, X*, gives 


Sx 2 c»x'2 


(5-28) 
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The operator L is thus a differential operator, of the second 
order with respect to x': 


L = -X' A + 1 AL 

bx 2 Sx * 2 


(5-29) 


Returning now to the old notations and adding the letter u to 
the velocity vector^ we get 


L = 


S 1 

u — + — 

Sx 2 Su^ 


(5-50) 


considering X*(t) as the velocity of point X(t) . The operator L, 
which is the subject of this example, appertains to the particular class 
of operators which make it possible to define statistical kinematics and 
consequently hydrodynamics, as will be proved. 

6. Random velocity. Hydrodynamic equations (refs. 5^ ^) • 

The problem involves the separation of what is position and what 
is velocity in the vectorial random function with six dimensions X(t) . 

In classical mechanics the velocity U(t) is the derivative of the 
position X(t) . In the present case, the assumption is made that the 
velocity is the stochastic derivative of the mean square of the position. 

By theory of random functions, it follows that, if the random function 
X(t) is completely or at least known locally, the function U(t) can 
be defined by operations comparable to those used in classical mechanics 
to deduce the velocity from the position. 

The preceding paragraph contained an example (equations (26) and (27)) 
of two random functions X(t) and X*(t) = U(t) linked by a relation 
of this kind, proving in a general way that there is an operator that 
satisfies the necessary relation: 


ARdu = 0 


(6-1) 
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ajad is such that the basic equation (5-6) (with original notations) 
teikes the form 

M + ^ Uj^M_ = aR (6-2) 

6t k 


If A is a differential operator, it can be expressed by an expan- 
sion in series (infinite or finite) in terms of the symbolical powers 

of the partial derivatives of which the coefficients are functions 

of X and u. The vectorial operator having for components occurs 

only in the first member^ by its scalar product with the velocity. 


The product G = RK can then also be decomposed in G = ppH (com- 
pare formulas (4-5) and (5-1)). The function H defines the correlations 
in the velocity field and function p is the mathematical representation 
of turbulent diffusion. 


Equation (6-2) must therefore play the fundamental part in the theory 
of turbulence. No attempt is made here at particularization; it simply 
is shown that it contains the general equations of hydrodynamics. 

First, the two members of (6-2) are integrated with respect to 
with due regard to (6-1). We introduce the density of the fluid 


p(x, t) 



t)du 


( 6 - 5 ) 


(or more accurately, a qmntity p normalized and numerically propor- 
tional to the density) and the relative mean of the velocity for a 
given position, the components of which axe 



( 6 - 4 ) 


We end with an equation of continuity 



( 6 - 5 ) 
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which simply expresses the fact that the velocity is the (stochastic) 
derivative of the position. Thus it is seen ttiat the velocity of the 
whole of the tiorbulent motion, defined qualitatively in paragraph 4, has 
the first of the qualitative properties of the hydrodynamic velocity in 
general. At the chosen turbulent scale (which is arbitrary) the mass 
is conserved. 

Now it will be seen that the velocity also satisfies the eqmtions 
of motion. Multiplying equation (6-2) by the velocity component 

and integrating with respect to u, gives for the first member 





Introducing the speed of fluctuation 


u'i = Ui - Ui 


and considering the equation of continuity, we put 


Tik = -Pu'i^’k 


( 6 - 6 ) 


The first member becomes 



the second member is written 



Uj^AR du = 


7, being a certain function of x which is deduced from the operator A 
and the probability density R. Hence the equations of motion 
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_8ui 

= ’'1 + p 




8t 


ik. 


8xk 


( 6 - 7 ) 


The play the part of the components of the density of the field 

of external forces (gravity, for example). The T^^^ are the components 
of the stress tensor (with fixed turbulent scale )5. They are, except for 
the factor -p, the components of the correlation tensors of the velocity 
components at a point. They can be measured directly by statistical 
methods, and the equations of hydrodynamics for a given turbulent scale 
can thus be verified. 

It should be noted that in this case the equations of motion of the 
ensemble (average) are involved. For the present, nothing about the 
behavior of the rate of fluctuation has been assumed. Later on the usual 
assumption will be made that the velocity satisfies the hydrodynamic 
equations and, more precisely, the Navier equations. It should be remem- 
bered that, based upon this hypothesis, Reynolds was able to establish 
equations similar to equations (6-7) for the mean turbulent motion. It 
is apparent that, without it being necessary to repeat Reynolds* calcu- 
lations, the statistical theory in question here is entirely different 
from that of Reynolds. It is probably more complete, but it still has 
not been pushed far enough to be verified by experiment, due to a lack 
of suitable hypothesis. 

As simple example of the fundamental equation (6-2), the case is 
chosen in which the function U(t) is itself differentiable (in mean 
squares, refs. I 5 , I 6 ). In this case, the statistical image point of 
the fluid has an instantaneous acceleration. If '^) represents 

the relevant mean of the acceleration, that is, of the velocity derivative, 
when the position and the velocity are fixed, it proves that the density 
of probability of X(t) and U(t) satisfies the equation 

^ i - 0 (6-8) 

Ot k=l OX^ k=l 

The operator A = - therefore linear, differential 

k=l 

and of the first order. The field of the external forces to which the 
fluid is subjected has, necessarily, for components, the quantities 


5rhis expression of stresses was originally given by Reynolds (On 
the dynamical of incompressible viscous fluids and the determination of 
the criterion. Phil, trans. Roy. Soc. CLXXXVI, part I, 123 ^ l895)^ 
proceeding from the Navier equations. However, the exact meaning of the 
Reynolds stresses is different from that of 


NACA TM 1577 


29 



(6-9) 


averages formed of the acceleration components when the position is fixed^ 
provided only that tends toward zero when the velocity increases 

to infinity. If^ in particular, the are not dependent on u, F^ is 

identical with and R satisfies the simple equation 



( 6 - 10 ) 


Once the 7^ are given, this equation defines R(x, u, t) from 
R(x, u, tQ^i. But the form of the linear transformation from R^x, u, tQ^ 
to R(x, u, t), that is, the probability of transition, depends upon the 
form of the original probability R(x, u, tQj, contrary to what happens, 
say for the equation of heat (type ( 5 - 19 )); this illustration does 
not appear to rest on hypotheses sufficiently inspired by reality to 
serve as basis of a turbulence theory. First of all, the starting point 
must be modified, as will be done in the following paragraph: 


7- Systems of molecules (refs. IJ and 14-): 

The statistical quantities to which the analysis of velocity records 
leads, represent only certain scales of turbulence, those which correspond 
to the ensemble of "vortices" whose dimensions are superior to a limit 
approximately fixed by the employed anemometer. Can the theory be changed 
so that all the possible scales can be represented simultaneously? It 
seems that it suffices for this purpose to start from the finest scale, 
that is, the molecular scale. The gas is therefore considered as a system 
of N molecules, and, to explain the method with as much simplicity as 
possible, the assumption, which is not verified for air, is made, that 
the molecules are identical, monatomic, comparable to material points 
subjected to central interactions. V(r) denotes the potential of the 
force of interaction of two molecules separated by the distance r. If 
r is great (with respect to the diameter of the molecules, which will 
not be introduced explicitly), V(r) is negligible. For the small values 
of r, V(r) expresses the repulsion of the molecules, generalized form 
of shocks. No other information about V(r) is needed beforehand, at 
least in a general theory. 

The motion of this system of molecules is controlled by the equations 
of dynamics. But the extreme complication of the trajectories of the 
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molecules prompted the replacement of rational mechanics by statistical 
mechanics. We prefer to introduce random mechanics where each molecule 
is a random point. By means of an ergodic hypothesis the random motion 
of a molecule can be considered as being statistically equivalent to the 
motion of the ensemble of the fluids in quasi-steady conditions, the 
ensemble of successive states of the visualized molecule replacing the 
ensemble of the simultaneous states of all the molecules . The molecule 
would therefore be the concrete image of the abstract random point serving 
up to now for representing the fluid. The molecular scale is a true ulti- 
mate scale of turbulence, separated, however, from the actual turbulent 
scales by a poorly defined but finite interval. While for experimental 
turbulence the concept of a random image point is a mathematical abstrac- 
tion, it is a natiiral idea and a starting point for the molecules. 


This idea is now explored but by a method slightly different from 
that discussed in the preceding paragraphs. The probability density 
fN(xi, • • of the positions and the velocities 

of N molecules simultaneously, rather than singly, is introduced. The 
gas appears then as a random point with three N dimensions, in a space 
of configuration, and no longer as a random point of ordinary and physi- 
cal space. For the time being, the notations x^_, X2, . . • ^ shall 

have a vectorial character and represent the system of the three coordi- 
nates of the molecules^ numbered 1, 2, . . ., N. 

The potential of interaction of the molecules of rank i and j 
is indicated by Vij, and the external force to which the molecule of 
rank , i is subjected, by 7^. If m is the mass of a molecule, the 
equations of motion of the molecules read 




dt 


4- 






(7-1) 


the summation applying to all values of j from 1 to N, when is 

assumed to be zero. Naturally V,- i = Vj,-. 

-L J J-L 

When these equations are written, the hypothesis is made that the 
velocity of each molecule is differentiable. It follows, that the random 

^In the preceding paragraphs the notation f(x, t) represented 
already the probability density of a vector x, of components Xj_, X2, x^. 

What is used only temporarily is the meaning of the subscripts, particular 
to this paragraph. A change in notation could be avoided only by compli- 
cations in writing which would be more harmful than useful. 
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point of 5N dimensions having for coordinates the ensemble of the coordi- 
nates of the molecules is represented by a random vector function of the 
configuration space of 5N dimensions doubly differentiable. This prop- 
erty which^ as stated before, is not necessarily true for the three- 
dimensional random point used previously, is the reason, according to the 
theory of statistical functions, why verifies the partial derivative 

equation 


at 1=1 ™ i=i 


(7-2) 


which expresses that fjj is ein integral of the eqmtions of motion. It 

has the same form as Liouville's equation of statistical mechanics. But 
it should he remembered that its original meaning is a little different. 
This equation replaces eind defines, in the space of 5N dimensions, the 
equation (6-2), and must now be exploited. 

The turbulence involves "particles" or "eddies" formed, at a given 
scale, by groups of s molecules, s being a very large number, but 
at the same time very small with respect to N. Let us Investigate what 
functional equation is satisfied by the probability density fg relative 

to the moleciiles of rank 1, 2, . . ., s. fg is defined by 




J 


f dx 


du 


N s-fl s+1 




(7-3) 


The fundamental equation (7-2) must be integrated with respect to 
^s-hl^ '^s-h1> • ' Xjj, Ujj. The operation is obvious, except for the 


terms in 


J 


; being function solely of u^^ and Xy the integra- 


tion gives a zero result if i and j are both superior to s . 

1 

If i and j are both inferior to s, then . 

® 6xj 

If J<s, i>s, the integration with respect to u^ shows that 
the result still is zero. 

If, finally, i £ s, j > s, we first can integrate with respect 
to xj^, Ujj, for k / 1. 
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The result is , where fg+1 represents the probability 

m ^ . Su . 

J i 

density relative to s + 1 molecules of rank 1,2,..., s, j. The 
last integration 



hf 


s+1 




dxs du^ 
J J 


1 b nbv^j 

mbu^J ^Xj 


dx^ du^ 
J J 


cannot be extended farther, because both 
of Xj. 




and fs+i functions 


Now it will be noted that in the s\m 



the terms are identical because the molecule j plays an anonymous 
part. 

The final result is 


at 


u. 


i=l 


afg 1 s 
— + - ^ 7i 


'ax 


m i=i 


a\ij^ 



■^s+l *^s+l ‘^’^s+l 


(7-^) 


This equation seems capable of serving as basis for a theory of 

tiirbulence. The scale there appears explicitly for the numbers i and 

s 

^ and the problem (not taken up here) consists in formulating a reasonable 

hypothesis which enables Is+q to be expressed with the aid of fg, so 
that (7-h) becomes a functional equation in fg. 
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For s = 1, (7-^) is, in a certain measure, comparable to (6-2) 

and reads 


Sfi 1 6fi N - 1 ^ rSV 


6t 


+ U-, — 


6x]_ m 


m <J Sxg 


^12 


•fg dXg du^ 


(7-5) 


Eqiiation (7-5) differs from (6-2) by the presence of the terms in 
fg* The presence of the second member relates (7-5) to (6-2), but its 

form makes it different from it. Equation (7-5) reduces to a functional 
equation in fj^ only if f2 is tied to f^^ by a suitable assumption. 

Such an assumption has been made by J. Yvon who demonstrated with its 
aid Boltzmann's fundamental equation on which the kinetic theory of gases 
is based. More recently it was tak,en up again by Born and Green, who 
also studied the case of s = 2, and applied it to the case of liquids on 
the basis of an assumption by Kirkwood binding f^ to f^. 

It seems that the case of large values of s has never been studied. 

Equation (7-5) like (6-2) and despite the not necessarily linear 
character of the second member with respect to fg, involves equations 

of hydrodynamics which are easy to write, if the method indicated in 
paragraph 6 is applied. Since these equations have so far not been used 
in concrete applications to turbulence, only these summary indications 
are given here. The application of equation (7-^) to turbulence raises 
difficulties which are far from being solved, or even stated but it is 
important to know that this equation exists. 


CHAPTER II 

CORRELATIONS AND SPECTRAL FUNCTIONS 


8. Introduction - correlations in space - homogeneity, isotropy: 

In this chapter we deal no longer with the general laws of proba- 
bility of turbulence . We limit ourselves to the study of the correla- 
tions of the law of the field H(u,u'j x,x',t,t'). Following are some 
preliminary remarks on the subject: 

Experiments fiurnish the correlations in form of temporal aversiges. 
Chapter I shows how it was possible to change their interpretation for 
converting them into stochastic averages . The assumption is made that 
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this operation is always possible. But for starting, there is no neces- 
sity for knowing whether it has been effected, since only properties of 
symmetry, of tensional character, are involved here. On the other hand, 
arriving at the dynamics of turbulence we shall see that the use of tem- 
poral averages leads to serious difficulties, which disappear when they 
are transfonned into stochastic averages, and our calculations will deal 
only with stochastic averages. Temporal averages are resorted to only 
when the calculations have reached the laws which must be compared with 
physical reality. 

The components of the speed fluctuations^ at point X]_,X2^xj are 

designated by Uj_, ^5- They are statistical quantities whose 

mean value is zero. Time plays no part at present and the parameter t 
in the formulas is disregarded. The correlations of the velocity between 
the two points x and x' = x + ^ are defined by a tensor having for 
components 


The turbulence is said to be homogeneous (within a certain field of 
space) when the are not separately dependent on the two points x,x', 

but only on their relative position, that is, on the vector | = x’ - x. 
Then Ra 3 (l) replaces Ra,3(x,x'). 

The turbulence is said to be isotropic at point x if it is homo- 
geneous in the vicinity of x and the tensor R^p is invariant to any 

rotation of the axes and any reflection. The form which isotropy imposes 
on the components of this tensor is discussed later. 

These two definitions concern only the second-rank tensor Ra 3 * 

They are "second-rank properties" of the statistical vector u^. They 
must be extended to certain third-rank tensors . Homogeneity and isotropy 
can be given a more complete and also more restrictive definition by 
extension to all possible moments of velocity components taken at any 
number of points. It is then more convenient to define the homogeneity 

^hey were called u'2_,u’2^u'^ in chapter I. But it is advantageous 
to use thereafter the notation u’ for other purposes. 



(8-1) 


The scalar is introduced also 



(8-2) 


a 
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as the invariance to translations of the law of probability of the 
ensemble of velocity vectors having for origin a certain number of 
arbitrarily distributed points^ and the isotropy, once homogeneity is 
achieved, as invariance to rotations and reflections . 

Without restrictive hypotheses the tensor has nine distinct 

components, fimctions of the seven variables X3^,X2,Xj,x'j^,x'2,x'^,t. 

If the turbulence is homogeneous, the nine functions remain, but of foxir 
variables only The isotropy is now written. 

To the tensor (5 ) we associate the scalar bilinear form 

y~ R .X 
^ — afi a p 

where and are two arbitrary vectors. This form is the mean 

value of the product of two scalar products ^ u^(x)X„ and 

^ Ujjj_(x' If it is invariant to rotations, it is an algebraic com- 
bimation, separately linear with respect to Xq. and Ya. of invariants of 
vectors Xq^,Yjj^,|q^ with respect to rotations. These invariants are 

y~xY y~iY = y~ ^ 2 

^ — a a a a z a a Z a 

There are therefore two scalars A(r), B(r), fionctions of r, so that 

H ''osVp ■ H Va + E ( z: (Z 5e’'p) (8-J) 

The identification shows that, if 6^^ represents the classical 
symbol of Kronecker, zero if p ^ a, equal to unity when p = a: 


^ap -^^aP ®^a^P 


( 8 - 1 .) 
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The tensor in the isotropy hypothesis, depends therefore solely 

on two distinct fimctions of two variables r and t. The notations are 
changed as usual and one puts with Karman and Howarth (reference 50): 


R „ = 
ap 0 


f - g 




(8-5) 


f and g being f inactions of r, t becoming unity for r = 0, and 
uq^ a simple function of t. 


If r = 0, R, 


aP 


is reduced to 




(8-6) 


It is seen that 


2 2 2 2 
V = ^2 = V " V 


U2U5 = u^u^_ 


\^2 = 


(8-7) 


Also, owing to the isotropy, the tensor R^^p is symmetrical. 
Lastly, it should be noted that the quadratic form ^ R X X 

OGp (X p 


is 


the first member of the equation of an ellipsoid of revolution. This 
shows that, if one takes, no matter how, a symmetrical table of numbers 
R^p, they are not, in general, components of a correlation tensor. These 

numbers must verify the inequalities which state that they are coefficients 
of the first member of the equation of an ellipsoid. In other words, the 
roots of the equation of the third degree 

S classical) must be positive. 


I^ap “ ^ (equation in 


A similar method permits the reducing of the components of the tensor 


’'oer “ + I) 


(8-8) 
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which will be needed later. This tensor is symmetrical with respect to 
a and p, and the trilinear form 




7 


is invariant to rotations. By applying the classical notations of Karm^ 
and Howarth, it is foimd, after a few calculations, that 


^aP7 ^ 




, 5 ( 3 - 1 / + — 5 I 


(8-9) 


a, and c being three functions of r and t. 

9. Properties of the functions f, a, b, c. Incompressibility^: 

Other fonns of symmetry less particular than isotropy can be vis- 
ualized^ as for instance, axial symmetry^ or invariance to rotations about 
a given axis, instead of about a point. The case of isotropy is that in 
which the depend on the smallest number of separate functions . 

The functions f and g are correlation coefficients. For example, 
f can be defined by taking two velocity components along the axis of 

(direction of velocity of the ensemble) at two points situated on a 
parallel to this axis (fig. 5 ) 


f(r) 


Ui^xi,xg,x^)ui(xi + r,xg,x^) 



(9-1) 


g is defined by taking two components still parallel to the axis of 
the Xj|_, but at two points located on a line perpendicular to this axis. 


^he authors (reference 50) ^se q, Z, k where we use a, b, c. 
But the letter k is also used in a just as classical although more 
recent fashion to designate the spectral frequency and will be used 
with this meaning in the present report. 
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For example 



(9-2) 


Experience indicates that near r = 0, the functions f and g are 
continuous and twice differentiable and allow a tangent to the ’’horizontal” 
origin (fig. 4). 

Theory confirms experience^ which unfortunately lacks precision when 
the distance r becomes small. Much greater accuracy is obtained by 
correlation measurements with difference in time (experiments by Favre, 
Report to the Intern. Cong, of Appl. Mech., London^ 19^8)^ because 

only one anemometer is used and the timing can be reduced as much as 
desired. 

Hence, one may write developments of the form 


f”(0) and g"(0) are negative quantities, possibly time functions, 
having the dimensions of the inverse of the square of a length. The 
length 


is called length of dissipation. 

The triple -correlation functions a,b,c have interpretations simi- 
lar to those of f and g. For example: 


f(r) = 1 + ^f"(0) + . . 




(9-5) 


g(r) = 1 + ^"(0) + . . 



(9-4) 
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It follows from formula (8-9) that, if has a well defined 

limit when | — >0, this limit can only be zero. Hence, a(0) ,b(0) ,c(0) 
are zero. In a general way, T^| 3 y is an odd function of 

Consequently, considered as function of r, c(r) is an odd function. 

Its development has no term in r^. Lastly, it is shown that it has no 
term in r, either. The coefficient of this term is the mean of 



u^_-^ = — lim ^ i — 

3 8xj_ 3 r — >0 r 


or 


i lim 
5 r — >0 


1 

r 



+ r 


l’^2^ 


13) - ui 5 (i 





But on account of the homogeneity extended to the averages of the 

third order, equal to u^5^xj^,X2,x^^ . Hence the 

limit is zero. Lastly, c(r) is an infinitesimal of the third order 
with respect to r. The same holds for h(r) and a(r). 


If the flow is incompressible, the various functions that charac- 
terize the correlations of isotropic turbulence are not independent. From 


the relation 



0 it is, in fact, immediately deduced that 



( 9 - 6 ) 


This equation is general. In the particular case of isotropic 
turbijlence it leads to Karman*s equation 


g 


2 8r 


( 9 - 7 ) 
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In the same manner, one then finds that: 


c = -2b a = -b - — ^ 

2 hr 


(9-8) 


So, of the five functions f,g,a,b,c which define the double and 
triple correlations of isotropic turbulence, only two remain distinct if 
the flow is incompressible. 


The incompressibility results, in particular, in the 
relation g"(0) = 2f"(0). Therefore, in terms of dissipation length A: 


f(r) . 1 . + 

g(r) - 1 - 4 + 


\ 




(9-9) 


Together with the length A, two other numerical parameters (time 
functions ) are frequently used to give an idea of the turbvilence . 

First, the correlation length 


L = 


f (r) dr 


( 9 - 10 ) 


is introduced; it depends upon all values of f(r) for 0 < r < <», while 
A depends only on the form of the function f(r) at the origin. It may 
be pointed out that L is expressed by means of the scalar 

R = u^(x)u^(x + l) = Uq (f + 2g). Owing to the incompressibility, 

f + 2g = 5f + r— = 2f + — (rf). If rf(r) approaches zero when r — >»: 
^r Sr 



2g) dr = 2L 


(9-11) 
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The intensity of turbulence in one direction, that of the eixis of 
for example, is the dimensionless quantity 



(9-12) 


where U is the mean velocity. The total intensity of turbulence is the 
ratio 


1 
5 

u 5 u 




u 


a 


yi(o) 




If the turbulence is isotropic, its intensity is the same in every 
direction. It is equal to the total intensity and to 

U 

In the case of isotropic turbulence the mean square values of the 
rotational components are directly associated with the quantities Uq 

and 7\. 


Introducing the components 




Su2_ 






Suj_ 

Sx2 


of the rotational, the averages CDp^, for homogeneous and 

and isotropic turbulence are computed: 



SU:r Sup 
2 — 2. _£ 
Sx^ 5x^ 
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We note that, for example: 


^7; 

— ^ = lim 
^^2 ^ 2 "^° 



^3 


(xi,xg,x^) 



= lim 
^2 — 


^2 


2u^(x3^,X2,x^)u^(x^,X2 + 


The first two bracketed terms are equal to Uq^ . 

P 

to Rjj(0,|2^0) = ^ 8> except for the third order 
Lastly 


The third is equal 



Likewise 


Sx2 y ^2 



2 V 

A2 
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Lastly: 




liin- 


^[^3(0^52. - 55) - Ri 2 ( 0 , 0 , i ^) - R25(o,52,oJ + Ri2(0»0,0’ 


Uq' 


,2 lim- 


12^3 


f - g. f 


2,7\n 


Finally 



2 2 2 2 
2 ^ + 2 ^ .+ ^ = 3 ^ 
7? a 2 ' a 2 a 2 


Natiirally , 




are eq\ial to 



10 . Spectral decomposition of the velocity: 

We know that the representative curve of a turbiilent velocity compo- 
nent u(t) in terms of time suggests the idea of an irregular periodic 
phenomenon. The analytical representation of u(t) is not a periodic 
Fourier series, but rather an almost periodic series or a Fourier integral. 

In the first case, u(t) is a sum of harmonics without common base 
period; by adopting the complex notation 


. . ^ io:^t 

u(t) = 

n=-cx> 


(10-1) 


the pulsations forming a succession of real numbers increasing 

with n. It can always be supposed that o)_n = a real 

quantity if the notation A^* designating the conjugate 


imaginary of A^. The series 


An 


must be convergent. 


MCA TM 1577 


kk 


In the second case : 


u(t) 



( 10 - 2 ) 


the function A(cd) being absolutely summable and such that A(-<d) = A*(o)). 

The numbers A^ or the function A(a)) depend upon the position of 
the anemometer. 

This representation is practical when the turbulence is steady in 
time, which enables a vigorous application of the ergodic principle and 
the calculation of the averages connected with u(t) in a time interval 
as great as desired. But, among the simplest problems involved, there 
is, first of all, that of a turbulence homogeneous in space (at least in 
a sufficiently restricted range) and developing in time. The pseudo- 
periodic character of the speed is then manifested in space rather than 
in time. (A detailed study follows.) 

The component u^ of the velocity of fluctuation at point X]^,X 2 ,x^ 

at a given instant can be developed in series or by Fourier integral, 
depending upon whether the spectrum is discrete or continuous 


u 


n^,VL2,n.^ 


^)+X2A2^ 

(t )e 


(10-3) 




(10-4.) 


ni, n 2 , and n^ are three integers that vary independently from 
-00 to +00 in formula (10-3), and on which the quantities Z. 


A 


(ni) ^(n^) 


a 


A, 


are dependent. In formula (l0-4), u^ is a 

triple integral extended over the entire space A of wave numbers (The 
wave number, inverse of a length, is the equivalent in space to the fre- 
quency, inverse of time. It is a vector), u^^ is a real quantity if 

Z.,(-Ai .-Ao.-Az.t'l = Z..*fAi ,Ao,Az,tV 
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These notations are of classical form. But, since the u are 

^ a 

steady, random functions of space, it is preferable to represent them as 
stochastic Fourier integrals 


In this formula, the hQ^(A,t) are random fimctions of 

with orthogonal increments (or noncorrelated) . Or, in other words, if 
A and A' are two different points in the space A, the mean of the pro- 
duct of the two increments dh^(A,t) and dhp(A',t) is zero. On account 

of the complex notations it naturally is a question of a product of 
hermitian symmetry. Hence 


If A' = A at the same point, this average is, in general, no longer 
zero. It is infinitely small of the order of dA. We put 


Now the spectrum of turbulence in the different cases (l0-5), (lO-^), 
and (10 -5) is defined. 


11. Spectral tensor and correlation tensor; 

Assuming spatial homogeneity it is now attempted to form, from the 
spectral decomposition of the velocity, the expression of the compo- 
nents abbreviated, of the correlation tensor 

of the velocity at a given instant. 

If the formiilas (l0-3) or (l0-4) are utilized, the spatial averages 
in a very great volume, physically limited but practically infinite with 
respect to microscopic turbulent lengths, must be used. It involves a 
cube with edges parallel to the axes, of arbitrarily large length 2a, of 
which the center, which can be placed anywhere becaiise of the homogeneity, 
is placed in the origin of the axes of the coordinates. 



dh^*(A,t) dhp(A',t) =0 if A' / A 


( 10 - 6 ) 


dh^*(A,t) dhp(A,t) = cPap(A,t) dA 


(10-7) 
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First, take the case of the Fo\irier series (10-5). To form 
involves, first, the product UQ^(x)up(x + ^), by associating an arbitrary 
term of the series which represents Uj^(x), or rather u^*(x), to any 
one term of the series that represents Up(x + |). To fonn the average 
in the volume V = 8a5, involves division by V and integration of x^ 

in space . The exponentials are restricted, and likewise their integrals 
in V, so that the q^uotients by V approach zero when a — ><». There 
is an exception for the terms of the series and u^(x + | ) which 

correspond to the same wave numbers, and which give the products 


2 





i 


e 





a 3 




(“p) 


p=l 


( 11 - 1 ) 


Independent of Xq^. They are therefore eq\ial to their mean value and 
consequently 


" u^^(x)up(x + I) = 

5 

1 ^ 

n^,n^,n^ 


( 11 - 2 ) 


This is none other than the Parseval formula for the almost periodic 
Fourier series of three variables. It naturally assumes conditions of 
uniform convergence which need not be defined here. 
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In the case of the Fourier Integral, u^(x)up(x + |) is a sextuple 
integral 


Ua(x)up(x + I) = 





5 3 

-i 21 Vp+i XI 

p=l p=l 


A' 


)(Xp+lp) 


dX dX' 


(11-5) 


extended over all values of Xj^, X2, Xj, 


X'2> and A'j. 


Integrating under the sign 


in the finite volume 


V 6uid dividing 


hy V, we find that, mathematically speaking, the mean value which is 
being sought is the limit for an infinite of 


i X. ^’p^p 

Z*(X)Z.(X')e ^ dX dX’ 

Y / a ' ' P' ' 

^ ■' AXA 


1 ^(X'p-Xp)'pcp 

P=1 


dx = 


n 1 > A' I 

/ ^ ^ ^ 

I / Z *(X)Zg(X')e ^ dX dX' 

^ c/axa 


sin ^X'j^ - Xj^j a 

^'1 - X 


sin 


(^’2 - >2) 


a sin ( X'^ - Xj 1 a 


X*2 - Xg 


X*5 - X5 


(11-4) 


Suppose Zj^(X) satisfies "Dirichlet's conditions" or, more gener- 
ally, is a limited variationj then the integral which figures in (U-4) 
has a finite limit which, when is continuous, has for value 
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This is Dirlchlet's theorem. 

The mean value of u^^(x)up(x + i) , the quotient of this integral by 

V = 8a5, thus approaches zero when a — > 00 ^ and the formalism of the 
Fourier integral does not furnish the spectriim. This conclusion is 
paradoxical . 

But every difficulty disappears when it is assumed that the turbu- 
lence might have only a line spectrum, incompatible with the so-called 
Fourier integral. But such an assumption does not seem reasonable. 

It seems more satisfactory to concede that the functions representing 
the tiirbulent velocity are too complicated for applying Dirlchlet's theorem 
and, more precisely, that their oscillations are too crowded to be repre- 
sented by functions with bounded variation. A mathematical process 
avoiding this difficulty is to follow. 

But an approximate argument can also be made . In reality , a is 
very great (with respect to turbulent wave lengths ) , but finite . Hence 
the mean value of u^j^(x)up(x + ^) is given by formula (11-4), without 

transition to the limit. However, a transition to the approximate limit 
can be made and assijmed that, with suitable accuracy 


UQ^(x)up(x + I) = ^ j Z^*(A)Zp(A)e ^ ^ dA (H-6) 


>^l] 


So, if a spectral function defined by 


1 > ^^p 




p=l 


dA (11-7) 


we get 




8«5 


Z^*(A)Zp(A) 


V 


( 11 - 8 ) 
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In a more accurate way, the cpQ^pCA) constitute the components of a 

tensor, the spectral tensor, transformed from FoLirier's correlation tensor. 
But the approximate formula (II-8) has the drawback of yielding a spectral 
tensor of a too restrictive form. In fact, there is no reason that the 
cp^p be the product of a function of subscript a by a function of sub- 


script 3, that is, the general product of the vector 



by 


itself. This new defect can be corrected by superposing a temporal average 
on the spatial average (II-8) . But it is much preferred to turn elsewhere 
and have recourse to the representation of the velocity by Fourier's sto- 
chastic integral ( 10 - 5 ). The product u^(x)up(x + |) is a double sto- 
chastic integral 


^a,(x)up(x + I) = 



^'pUp+Sp) 


dh^*(A)dhp(A') 


( 11 - 9 ) 


On averaging, it is seen that the qviantity dh *(A)dh (A') is 

QC p 

involved which is, as explained earlier, zero when A' / A, and of the 
form cpo,|3(A)dA if A' = A. Hence simply 




_1 


(11-10) 


The correlation tensor ^0,13(1) is the Fourier transform of the 
spectral tensor <Pcc| 3(A). The here have the general ty desirable. 

They are subject to the two following conditions : 


(a) 9g^p(^) has hermitian symmetry 


V = 'P3a* 
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(b) Whatever the complex numbers the hermitian form 

V|3*V 

a, 3=1 


which is real according to the first condition^ cannot be negative. 

Any system of numbers cp^p constituting a tensor and satisfying 
these two conditions can, a priori, be used as spectral tensor. 


12. Spectral tensor of isotropic, incompressible turbulence: 

I. Suppose, with Heisenberg (reference 25 ), that the turbulence is 
isotropic. Then, as supplementary condition, the flow is incompressible. 

The method used to express the isotropy of the correlation tensor 
is applied to all tensors, whatever its meaning. The spatial tensor is 
isotropic when two functions A(k), B(k) exist depending, moreover, on 
t such that 


cPap(^) = A(k)A^Ap + B(k)6^3 


( 12 - 1 ) 


2 2 2 2 
k denoting the quantity + >^2 + . 

If the flow is incompressible, the functions A and B are formed 
by a relation equivalent to Kerman’s relation (9-7) between the functions 

of the correlation f and g. We write, in fact, that ^ ^ = 0. 
Since 




= i 


P=1 


A dh 


a 




a a 


(12-2) 
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we must have identically 


Z A dh = 0 

a a 


(12-3) 


hence, with 3 being fixed 





= 0 


When this condition is applied to the components (l2-l) of the 
isotropic spectral tensor, it is seen that 

Ak^ + B = 0 


We put 


hence 


A 


_F 



m L 

V = ^(^03 



(12-4) 


F(k) is the spectral function introduced by Heisenberg. It will be 
interpreted later. 

II. With Kampe de Feriet (reference 26) we first express incompressi- 
bility without disturbing the isotropy. 




are subject to the condition that the hermitian form 


♦ • Z V3* foB 


(12-5) 
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is positive or zero. It is positive, except when the vector is 

proportional to vector since 


on accoimt of the incompressibility condition (12-5) • 

The hermitian form \|r is now reduced. According to classical theory 
and the equation in S, three real numbers S]_,S2^Sj and three linear 

forms 


the three vectors of components two by two orthogonal. 

If t — 0^ the numbers Sj_,S2^Sj are positive or zero. 

But, there exists a particular vector 7^, transformed from vector 

by (12 -7)^ so that \lr = 0. This can happen only when the three terms of 
which \|r is the sum, are zero at the same time. Hence two separate 
hypotheses 



(12-6) 




can be found with which \|r can be expressed in the form 



(a) S2 — 0 — 0 






( 12 - 8 ) 
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being the function of A. This form of cp^^, similar to (11-8), is 
too restrictive; 

(b) S3.0 

♦ = % S 2 (12-9) 

We write rather 


♦ ' ®1 H + S 2 ^ 

oc3 a$ 

and identify with the original form (12-5) of f. It is seen that 


V = Sia^*ap -h S2b^*bp 


( 12 - 10 ) 


This is the most general form of a spectral tensor in an incompress 
ible medium. In more simple form 




= b 


! 

a 


There are two vectors forming with vector a trirec 

tangular trihedral and such that 


'Pap = + ^'a"^'p 


( 12 - 11 ) 
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By Pythagorean theorem 




r.'^nr 

a a 




(12-12) 


with 


a'2 




,.2 



, 2 
3 


of 


By this formula, the b'^ can be eliminated from the expression 


'Pap ■ 



(12-15) 


Finally we put 


^,.2 ^ F(k) 
4«k^ 



This is the canonical form which Kampe de Fe'riet has given for the 
spectral tensor of an incompressible fluid 




ap 



( 12 - 14 ) 


with 


Z =a(>')V ■ ° 
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The isotropy stipulates that = 0 , so that formula (12-4) is 
obtained again. 


15. Energy interpretation of the spectral function F(k): 

Except when stated otherwise, the isotropic turbulence is assimied 
to be in the incompressible state. 

Hence 

2 2 2 2 
ui^ = = uj*^ = Uq'" 

The total energy of the turbulent fluctuations per unit mass is 


E = 



1 

2 







(15-1) 


But 


Hence 


9 (4) 

Z — ^oa 


= 

2nk^ 


E = ^ r Eilil dA 

Ja k2 


(13-2) 


(13-5) 


In passing to the polar coordinates in the space of the wave numbers 
one finds that 


E 


loo 


F(k)dk 


(13-4) 
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Which, after minor changes, gives the amount of the energy of the 
turbulent fluctixations per unit mass in the "sphere of the wave numbers" 


E,, = 


F(k' )dk' 


(15-5) 


k’ being an integration variable. 

The ’’small values of k" correspond, of course, to the large scales 
(large vortices). 

The next problem is to find the expression of the energy e dissi- 
pated as heat per unit mass. While the preceding calculations were by 
nature strictly kinematic, the dissipative function here must be so intro- 
duced that it yields equations of motion, that is, Navier equations; 
hence, an assumption associated with the dynamics of turbulence, which 
is discussed in chapter III. 

With V denoting the coefficient of kinematic viscosity, the energy 
dissipated in heat by the molecular motion is 


V 



+ 




(13-6) 


This is to be expressed in spectral terms and averaged. 


we get 



If 


NACA TM 1577 


57 


and 



which, squared, reads 







dh^*(A) + 


A, 


a 


dhp*(A) 


A' 


+ ^'a (^5-7) 


To obtain the averages, 


involves operation under the 



The 


result is zero, except when the points A and A' coincide in the 
space A of the wave numbers . It leaves 




An 


dh, 


a 


''a''3 




+ A„A« |dh„* dhp + dhp* dh, 


^a] + ^a^KI^} 


+ \t^3('Pa3 + ^3a) V^'P33]‘^^ 


(13-8) 


Consequently, if e is the mean dissipation of energy per unit mass 


e 



9, 


oux 



V^ 3 ^a 3 


dA 


(13-9) 
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This formula supposes neither incompressibility nor isotropy. If, 


now, the turbulence is incompressible, ^ ^)wing to the 

ap 


isotropy, YZ ^oo. " ‘ 


2itk 


Hence, by integrating in polar coordinates in space A 


e 


2v I k%(k)dk 


(15-10) 


e being a physical quantity, hence finite, it is seen that the 
preceding integral is finite, which gives a first limitation of the 
possible forms of the functions F(k). 


Ik. Relations between spectral function F(k) and correlation func- 
tions f(r) and g(r): 

The spectral tensor was defined by the formulas 







which assume homogeneity only. Incompressibility and isotropy are to be 
added. 

According to section 8 






0 



+ S6o,p 



8r 


and according to section 12 




F(k) / ^a^p\ 

4.k2 V-P ■ k2 ) 


NACA TM 1577 


59 


By Fo\irier's reciprocity formvilas 


5 

n ^ 


the integral being extended oVer the entire space. The general formulas 
for defining F in terms of f or g, and f, g in terms of F must be 
part i c ular i zed. 


It is pointed out that 






OJOL 


is the scalar invariant of the 


spectral tensor 


cp^p ("trace of tensor" or contracted tensor). Hence it 


seems logical to compare with F the analogous invariant of the correla- 
tion vector defined by 


R(r) = ^ R^ca = + 2 g) = UQ 2 ( 3 f + rf) 


Thus the following problems must be solved: 

(a) express F(k) by means of R(r)^ and conversely 

(b) express f(r) and g(r) by means of R(r) 

(c) express f(r) and g(r) by means of F(k) 

Also to be defined in precise manner is what is called the trans- 
versal (lateral) and the longitudinal spectrum of turbulence. 

(a) Relation between F(k) and R(r): 




Vp 



^ I 
p p 


F(k) 

2jtk^ 


dA 


(lk-2) 


R(r) = 


(a^ 


dA = 


e 
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To reduce the triple integral to a simple integral, simply select 

a new axis of the perpendicular to the plane ^ ^^p “ 

resiilts in the Integral 


R(r) = fe^’^^5 dA = — f ®F(k)sin 6 dk d6 dcp 

Ja 2rtk2 2n 


(14-3) 


Which, after a short calculation, yields 


^00 


R(r) = 2 


Jo 


sin rJ 
rk 


^(k) 


dk 


(14-4) 


Conversely 


= — r sin rkR(r)dr 

k « Jo 

(h) Relations between f(r), g(r), and R(r): 

The integration of the differential equation yields 


(14-5) 


rf ' 


+ 5f = ^ 


(14-6) 


or 


^(rM = ^ 


B 
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with the condition f(0) =0. It is 



— I r'%(r')dr' 

r5 Jo 


(14-7) 


2 

R(r) having for limit -Juq , when r tends toward zero, it con- 
firms that the integral is equivalent to , and that f(r) tends 

rather toward unity. 

From that, it is deduced, by formula (9-7) that 


UQ^g(r) = |R(r) 


1 

2r5 



%(r ' )dr ' 


(14-8) 


(c) Relations between f(r), g(r), and F(k): 

A simple calculation gives 


UQ"^f(r) = 2 F(k) 

b V r5k5 


sin rk cos r 


k\ 


A2 ! 


dk 


p / , I T./, \ /cos rk sin rk sin rk\ 
UQ^g(r) = F(k) — — + dk 

Jo \ r^k^ r5k5 rk / 


(14-9) 


From that the first terms of the limited developments of R(r), f(r), 
and g(r) for small values of r are derived in terms of the dissipated 
energy e 
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R(r) = 5 Uq 2 1 - 

\ VUq'^ 


f(r) = 1 - — r^- 


1_2 e 


50 VUq^ 


(lU-10) 


g(r) = 1 - ^r 2 _§_ + 

1*5 2 
^ VUq 


Comparison with (9-9) indicates that the dissipation length A is 
related to the dissipated energy e, the viscosity v and the kinetic 

■5 p 

energy of the tirrbulent fluctuations E = throiigh the formula 


F(k)dk 


a 2 = lOv^ = 5- 


E 


k^F(k)dk 


(14-11) 


15 . Lateral and longitudinal spectrum: 

The spectral measurements do not provide the spectral function F 
directly, but merely the spectral terms corresponding to certain simple 
associations of velocity components, for the observers placed in the 
particular relative positions . It concerns longitudinal components of 
the velocity (parallel to the velocity of the main flow), at points 
placed either on a parallel to the axis or on a perpendicular to the 

axis. The corresponding correlations are UQ^f(r) and UQ^g(r). 
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According to the original notations of G. I. Taylor 





\ 


COS 05rA(a))d40 


> 




cos oirB(cD)d4D 
y 


where A(o)) and B(a>) are termed the longitudinal and lateral 
functions. They now must be expressed by means of F(k). 

Since g = f -h — f ’ ^ elementary calculation shows that 


B(o)) = “[a(o)) - o)A^(a3)J 


and at the same time that 


R(r) = 2 



rk 


OCX) 


Jo 


cos osr 


[A(oi) + 2 B(o))] da) 


Elementary Fourier transformations yield 


A(cd) + 2B(cd) 



A and B are determined by the two equations 


2B = A - <juA' 


Jo k 


(15-1) 


spectral 


(15-2) 


(15-5) 


(15-4) 


A + 2B = 2 


(15-5) 
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The result is 




(15-6) 


CHAPTER III 


DYNAMICS OF TURBULENCE 


16 . Introduction: 

To construct a dynamics of turbulence, that is, to set up the dif- 
ferential or finite laws which govern the development of the statistical 
quantities characterizing the turbiilence in time and space, it is neces- 
sary to start from elementary laws and apply the statistical methods to 
them. The most natural idea, the only one which actually produces con- 
crete resiilts, consists in utilizing the Navier equations. To what 
extent are they applicable to turoulence? The turbulent motion is always 
a macroscopic motion with respect to a finer scale motion, and, at the 
limit, with respect to the molecular disturbance. Therefore it is rea- 
sonable to believe that it satisfies the equations of the mechanics of 
fluids. The next step is to find the solution of the Navier equations 
which, for certain limiting conditions, have the turbulent aspect, and 
calculate the particular averages from these solutions. Unfortunately, 
it is rather difficult to define these limiting conditions. So, the 
remaining resource is to examine whether , among all the possible solutions 
of the Navier equations, there exist any sufficiently complicated for rep- 
resenting the turbulence, without attempting to determine them logically 
by the limiting conditions. But then a new difficulty arises. Solutions 
for the Navier equations are known only for simple conditions which are 
far from resembling turbulence. In other words, while conceding their 
validity, we practically do not know how to solve them. 

Since it is not acceptable to take the averages on the solutions of 
the Navier equations, it is attempted to take the averages on the differ- 
ential equations themselves and to write the differential equations which 
verify the statistical quantities. Since this method tos produced results, 
it is set forth in the simplest case, that of the Karman correlation 


tensor 
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17 . Fundamental equation of turbulent dynamics: 

It is assumed that the total velocity U is uniform and along the 
axis of The Navler equations read then 


^ i: = .i ^ ^ V Au, 

at dxi 3 “^axp p ax^ 


(17-1) 


in the absence of external forces, p is the density, usually assumed 
constant, and p the pressure. These equations must be supplemented 
by the equation of continxiity 



(17-2) 


The notations are abbreviated by representing the velocityT at 
point x' = x„ + 5 by u' = u„(x' ) = u^(x + § ). Multiplication of 

Qw CX QC CL CL (X . 

(I 7 -I) by u* and summation with respect to a, gives 

CL 


^ u' + U / u' + ') u' Uq^^ = 


a "a 




, ^ 


u' — + V 


a^u 

u' Si 

op” 


(17-5) 


hn chapter I, section 6 , represented the velocity compo- 

nents of the whole and u' 2 ^u' 2 ,u'^, those of the velocity of fluctuations. 

The over-all velocity having then U,o,o, as components, the notations are 
changed; Uj^ represents the velocity fluctuations at point x, and u'i 

the velocity fluctuations at point x ' . 
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Permutation of the points x and x' yields a similar equation, 
which is added to the preceding one 


6 

6t 


^ a a V \ V V *^^^3 ^^ 3 / 






a/ a3 


u 




u' 


6xp^ Sx'p^y 


( 17 -lf) 


Taking the averages of the two members it is then assumed that the 
turbulence is homogeneous, nothing more. The scalar product ^ ^a^'a 

is a function of 1^^, and is not individually dependent on. x^ and x'^^. 
Consequently 




Sx 




6x' 




The term containing the over -all velocity U disappears. The 
following term can be written 




6u 


U' U_Uo + / U u' U'd - N 

o3 o3 ' 3 3 ' ^^3 ^ 


a a 


6u' 


3 • p a 


^ ' u_u' 

a a 
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The last two terms disappear as a result of the equation of continuity. 


Y1 is written :8 

OXp 


^ u' U^Uo = 
— a a p 

P a 


P din a 


Likewise, as a result of the homogeneity, ^ 


u„ u' u'o becomes 

°-A-y ' P 

ox p 




Since ~ system of these two terms is equal to 


51 7 ^ 5lT„»a(0 


P ^Ir a 


(17-5) 


The viscosity term reads 


2v ^ ^ = 2v A ^ - 2v AR (I 7 - 6 ) 

P a a 

where A is the Laplacian symbol in three-dimensional space. 

Lastly, it is shown that the pressure terms disappear. This is done 
by extending the definition of the homogeneity to the averages containing 
the pressure, with due regard to the Incompressibility. The averages such 

as pu'(x functions separately of x and x* but merely of the 

difference | = x' - x, so that, on these fimctlons, 

ox ' ox 

Q- 

See the definition and the properties of the tensors R^^p 
'^aP7 chapter II, sections 8 and 9- 


and 
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Therefore, u'^ being function of x' 




u' 


^ = 


a- 


Bx, 


E 


a 


a 


= E 


a 


Su* 

Sx' 


= 0 


a 


(17-7) 


Likewise for the other pressure term. The final resvilt is the 
fiondamental equation (reference 50 ) 


|:EV = 2vAXIV+25lTpE’'a3a (^T-8) 

Ot a, a 3 O' 

It was written with the aid. of the components of the tensors 
and 1 a,p 7 * Introducing the scalar R = ^ ^rrrt.^ vector T which 

has for components 


^3 ^ ^a3ct 

a 

The equation (17-8) assumes the form (reference 2 ) 


(17-9) 


^ = 2v AR + div T (l7-10) 

Bt 


l 8 . Case of isotropic turbulence: 

The equation ( 17 -IO) assumes Incompressibility and homogeneity but 
no isotropy. If the turbxilence is isotropic, it is known (l4-6) that 


f being the first of the two Karm^ correlation functions. 
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According to the properties of (section 9) 


’'S ' + 2a) = Uo\(c 


( 18 - 1 ) 


The divergence of this vector is equal to 


div "t" = ^ 1 

Sr 


rJfc' + !£' 


(18-2) 


The fundamental actuation (17-IO) is now transcribed: 


r2 Sr 



2v S /^Sr\ 
j,2 Sr V Sr/ 


r2 dr 




In fact, it is remembered that the Laplacian of a fimction cp(r) in 
n-dimensional space is 


Ancp(r) = cp" + ^ ~ ' 


(18-5) 


Multiplication by r^ is followed by integration with respect to r. 
It is easily verified that all the terms cancel out for r = 0. Bearing 
in mind that 




one arrives at the equation 



(18-4) 


(18-5) 
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This is the Karm^-Howarth eq.uation (ref. 30 ) in its classical form. 

The general equation (17-IO) has the form of the equation of heat 
propagation in ordinary three-dimensional space, the quantity being 
propagated is the scalar R, the invariant of the double correlation 

tensor. The "second member" div T is tied to the triple correlations 
as R is to the double correlations . The result is an equation of 

partial derivatives between two \mknown functions R and div T. For 
reasons which appear later on, R is usually considered as the principal 
unknown. The problem of solving (17-IO) consists in defining R by means 
of suitable physical hypotheses on the triple correlations . The equation 
with two unknowns finally arrived at is the result of the nonlinear char- 
acter of the NavLer equations . Its mode of operation will be explained 
later. 

The "isotropic" equation (18-5) leaves no trace of the original 
properties of tensorial symmetry. Its purpose is to connect the func- 
tions f and c, which can be measured directly. If it is borne in mind 
that 


f"(r) + ^f’(r) = A5f(r) 

> 


c' 


r hr 



(18-6) 


then one may write: 

(uo^f) = 2v A5 (uo%) . ^ |:(rV) (18-7 ) 


It is an "eq\iation of heat," isotropic naturally, but in a fictitious 
five -dimensional space. This remark is interesting, but it seems arti- 
ficial and without real physical signif icaince . 

Nevertheless, an un\isual property of the functions UQ(t) and f(r,t) 

can be derived from it. Multiplying by r^dr followed by integration 
from 0 to oo yields 
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Hoo 

r^f dr = 0 

° Jo 


(18-8) 


II ^ ll 

provided that r f and r^c tend toward zero when r — > 00 ^ which seems 
reasonable, and difficult to check. Consequently 


'^00 

UQ^(t) r^f(r,t)dx 

Jo 


(18-9) 


is a nvimericaJ. constant, independent of t, d\iring the development of the 
tnrbulence. It is the Loitsiansky invariant, which connects the energy 

of fluctuation E = "to "the correlations f. 

It also should he noted that, since the second term of ( 17-10) is 

the divergence of 2v grad R + T, its integral, on a sphere of center 0 

ajid infinitely large radius, is zero. Consequently, the integral of 

8R , , 

— xs zero, and 



dxi dx2 dx^ = 


(18-10) 


the triple integral being extended over the entire space. If the turbu- 
lence is isotropic, it leaves 


r2R(r,t)dr = 


(18-11) 


But, according to the expression of 
beginning of this paragraph: 


R(r), remembered at the 


r^R dr = Uq2 




0 r 


lim 


(rJf) 


72 


NACA TM 1377 


If 


r^f ' 


>•0, r^f tends tcwaxd zero also, and this limit is zero 


poo 

t2r dr = 0 


u 


0 


(18-12) 


p 

This proves that the fimction R, which is equal to 3^q for 

r = 0, and to zero for r = », takes negative values for sufficiently 
high values of r. If formiila (l8-12) is approximated to formula (lA-5) 
which expresses F(h) by means of R(r), it will be found that, pro- 
viding R(r) tends sufficiently fast toward zero at infinity, the 
development of F(k) in powers of k starts with a term in k^; near 
k = 0 (large size eddies) 


F(k) = Ck^ 


(18-15) 


with 


= r^(r)dr 


6it 


The constant C can be expressed by f, because, according to (lA-7) 


Uq^ I r^f(r)dr = / r dr / r'%(r’)dr' 

Jo Jo Jo 


In order to change the order of integration, the integrals of the 
second member are written in the form (fig. 5) 


lim I r" dr" 1 r'%(r')dr' = lim I r'%(r')dr' / r dr 


r — ^ 00 


T > 00 


i lim r^ / r'^R(r')dr' 

2r-^~ Jo 


1 r''4^(r')dr’ 

2 Jo 
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Now: 



r '%(r' )dr' 


u^T^t{r) 


So, if f(r) 
not only r^f but 


tends toward zero at infinity fast enough so that 
even r^f ' approaches zero^ then 


2 poo 

C = I r^f dr (18-14) 

5^* Jo 


To a factor — , C is therefore identical with Loitsiansky ' s 
5« 

invariant ( l8 -9 ) • 

These resiilts could be extended to include the nonisotropic 
turbiilence . 

Lastly it should be noted (ref. 2) that, in the case of isotropic 

turbulence, the vector T defined by (l8-l) is radial, that is to say, 
that the Tp are proportional to the |p. In effect: 



Uo^(c + 2a) 


(18-15) 


Hence one may put: 

~T^ = T grad r (l8-l6) 

T(r) = Uq^(c + 2a) being a scalar. 

19. LocaJ. form of the fundamental, eqiiation; 

The problem is to ascertain what the equation (17-IO) becomes when 
r — >0. 


74 


NACA TM 1577 


2 ] 

It is known that f(r) — >1. As f = 1 - — — + . . ., f" + — f 

2X2 

5 

approaches — 

X2 

Lastly, as c is of the third order with respect to r, c' + — 

r 

tends toward zero. 

The final equation, due to G. I. Taylor, reads 



Ur 


-lOv- 




(19-1) 


Recalling that, according to paragraph l4: 



x2 = 10v£ 
e 


The formula (I9-I) becomes therefore 


9 



(19-2) 


It establishes an elementary relationship between the dissipation 
of energy e by viscosity and the decay as a function of time of the 
energy of fluctuation E; the turbulent energy of fluctuation is totally 
dissipated as heat by viscosity, this dissipation following a process to 
be analyzed in chapter V. 

The foregoing demonstration stipulates isotropic turbiilence . But 
the resxilt (19-2) is valid in more general cases. In fact, (19-2) can be 
derived from (I7-IO). First it is known that R(0) = 2E. The connection 
between A R and e is established by the intermediary of the spectral 
fimctions cp q(X). If the flow is incompressible, formula (l3-9) reads 

Oup 


being function of the sole variable t, ^ can be replaced by 

ot 

— , a notation which represents an ordinary derivative with respect to 
dt' 

time (not a derivative with respect to motion which would be without sense 
here ) . 
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e 


V jfk2 ^ cp^(\)d-K 


(19-5) 


but on the other hand^ 


R(l) = XI WO 



q>^(-K)d-\ 


consequently : 


AR = 


k^e 


p=l 





(19-4) 


and in particular, when W approaches zero, v/^ is reduced to 



To prove that div T = 0, at the limit, the simplest way is to revert 
to the Navier equations. 

It is plain that equation (19-2), limit of (17-IO) when r — >0 can 
he obtained without going by way of the finite values of r. Simply mul- 
tiply the Navier equations by U2,U2,u^, add them up and average. Assinning 

the disappearance of the pressure terms to have been proved, it is 
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The term — — u ^ is written 
2 dt ^ 


dE 

dt 


The second member. 


according to the foregoing proof, is equal to -e. Hence it must be 
shown that 'S u u„^^ = 0. This q\iantity reads 

^ O' 


1 

2 



1 o 2 

— / ^ ^ Uq 

O ^ CL P 


-1 Hu 

2 a ' 


a 


Sup 

Sxc 


(19-6) 


By virtue of the homogeneity, u^^Up is not dependent on x, and 
the first term is zero. 

Owing to the incompressibility, the second term is zero. The first 
member is therefore zero, which provides the proof. 


20. Solution of the fundamental equation, when the triple correlations 
are disregarded: 

When the triple correlations are discounted, equation (I7-IO) becomes 
the equation of heat 


^ = 2v (2O-I) 

at 


This hypothesis, which mathematically is convenient, is physically 
quite difficult to jiastify. Its accuracy increases as the viscous 

effects become greater (prevalence of term 2v AR over the term div T ) , 
or as the Reynolds numbers decrease (whatever their definition). It also 
wovild be verified if the velocity components followed the Laplace- 
Gaussian law. And, what is more interesting, it may be added that it 
corresponds to a state of turbulence in which the forces of inertia are 
negligible against the forces of viscosity. 

Equation (20-1) must thus be resolved knowing R for t = 0, and 
being aware that, if r— >0, R has for limit a finite quantity Juq (t) 
temporarily considered as known. 
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In isotropic turbulence, R depends only on r, t, and (20-1) reads 

( 20 - 2 ) 


dR „ 2 

^ = 2v + - 

at y^j.2 r ary 


To satisfy the limiting conditions, elementary solutions of the form 
R(r,t) = Rj^(t)Rg(s) s = ^ (20-5) 


\J 8vt 


are necessary. Hence 


1 — ( 

1 

-p 
1 — 1 

R"2(s) 

R^(t) 2R2 (s) 

2 


+ (s + -)R' 2 (s) 


(20-3) 


Q 


The two terms of this equation have a constant value - There- 
fore, R 2 satisfies the linear differential equation 


sR^ 


{"2 + 2p + i)r' 2 + QsRg = 0 
or, reduced to classical form, by the transformation, 


(20-A) 


sRg = e"® H 


(20-5) 

Quantity H confirms in fact the Hermitian differential equation 


H" - 2sH' + (Q - 4)H = 0 


( 20 - 6 ) 
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Since R(r,t) must cancel out for r infinite, hence also H(s) 
for s infinite, Q - 4 is necessarily an even integer 2n. In that 
case H is the nth Hermitian polynomial: 


Hn(s) = (20-7) 

ds 


and 


R2(s) = (20-8) 


But R 2 must have a finite limit when s — K), which requires that 
H^(s) he an odd polynomial, hence that n be an odd integer: n = 2p - 1, 

Hence the ele- 


Q = Ap + 2. The function R 3 ^(t) is then equal to 
mentary solution: 

R(r,t) = 





(20-9) 


The general solution is a superposition of elementary solutions : 

the constants Ap being so chosen that the series converges and R(0,t) 
is equal to 5 Uq^. 

Examples: 1. Limited to the term in p = 2 it is found that, by 

interpreting the Hermitian polynomial, Hj(s), and denoting an inde- 
terminate constajit by A: 
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t2 

hence 

_ r2 

f(r,t) = e ^ 

This example is rational. It is compatible with the hypotheses 

made at that time. In particular^ r^ approaches zero when r — ^ 00 ^ 
whatever the non-negative number m may be. However, f = 0 for 
t = 0, except when r = 0. At the initial instant the turbulence is 
concentrated at one point from which it ultimately spreads throughout 
the entire fluid. Such a turbulent structure seems rather difficult 
to conceive when assimiing it to start from t = 0. An interpretation 
will be given later (section ^k) in connection with the final phase of 
decay of turbulence behind a grid. 

For p = 1, one would have 


_5 

UQ2(t) = At 2 


( 20 - 11 ) 


r^ 

E(r,t) = ^'8vt 
t2 


f 





u 


0 


( 20 - 12 ) 


f wOTild not converge toward zero when r — which proves that this 
example cannot be suitable for a real motion. 

2. (K^rm^-Howarth (ref. 50)- Suppose the correlation function f 
does not depend separately on r and t but only on the vari- 
able s = - 

\J 8vt 


Since 
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it then yields 


R = UQ^(t)j5f(s) + sf'(s^ 


(20-15) 


so that R is the product of a function of t by a function of s. 
Thus R is an elementary solution of the problem treated in the first 
example, corresponding to a chosen Hermitian function 

The function f(s) satisfies a differential equation similar to 
(20-4) written by K^m^ and Howarth. But a great advantage accrues 
from the use of function R and reduction to Hermitian polynomials . 

Remark. According to (20-10), R(r,t) is a sum of "elementary 
solutions" corresponding to various values of the integer p, starting 
from p = 2. For each one of these solutions, except for p = 2, it 
is verified that the invariant of Loitsiansky (l8-10) is zero. For 
p = 2 it has a finite value . 


21. Solutions involving a similarity hypothesis (see ref. j): 

If the triple correlations are no longer negligible, the fundamental 
equation cannot be solved completely. It is assimied that the correlation 
functions f(r,t) and c(r,t) are not separately dependent on the two 

variables r,t, but solely on \lr = ^ where l(t) is a length in 

Z ( t ) 

terms of time. That is to say, that at each instant t the correlation 
curves are superposable , by means of a unit change on the axis of r. 
However, it is recognized that this hypothesis, called "total similarity" 
is often too specific, and it is therefore replaced by a "partial 
similarity" which is verified only in a finite interval Zp < < ip • 

For total similarity, the Loitsiansky invariant reads 




I '00 1 

2(t) r^f(r,t)dr = Uf.^(t)z5(t) 'Jr^f(t)dt 


Jo 


( 21 - 1 ) 


Consequently, UQ^(t)Z^(t) is a constant independent of t in the 

course of the motion. This is no longer true in the case of partial 
similarity. 
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To transform the fundamental equation, it is possible to put 


R(r,t) = 5Uo^(t)a(>|/) 


T(r,t) = 3Uo^(t)p(t) 


( 21 - 2 ) 




UqZ 


a(i|r) is a fvinction that takes the value unity for >|f = 0; p(\|r) is 
a function that cancels out with \|r, is a Reynolds number in terms 

of the time, and T is the scalar which, according to (l8-l4), defines 
the triple correlations. 



Considering the relation (19-I) which affords 

dt 

and of the length A, the fimdamental equation 


in terms of 


— = 2v /SR + div T 
8t 


is easily transformed to 


^(z^ + a"(ir) + ^'(t) + fp'('l') + 

4v ' ' ^2 t 2—1' L i|/ 


= 0 


(21-3) 


or, schematically arranged: 


U'l('t)3i('l^ ) + ci2('t)32('l^) Q,i|.(^)3l4.(’^) - 0 (21-4) 


Equation (21-3) is a sum of foiir terms each of which is the product 
of a function of t by a function of ilr. How can an equation such as 
(21-4) be resolved? 
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The Oj_(t) are considered as the coordinates of a point A in 

four -dimensional space ^ this point describing a curve (a) in terms of 
the parameter t. The same applies to the It results in two 

curves (fig. 6) such that the straight lines joining the origin to any 
point A of the first curve and to any point B of the second are 
always perpendicLilar . These curves must fit two orthogonal complemen- 
tary subspaces in foirr -dimensional space. Hence there are^ a priori^ 
three possible cases: 

1. (a) is a line passing through the origin. (b) is in the com- 
plementary three-dimensional space orthogonal to this line. 

2. (a) and (b) are in two completely orthogonal planes passing 
through the origin. 

5. (a) is in a three-dimensional linear space, (b) is the orthog- 
onal line. 

In the first case, the aj^(t) are proportional to constant num- 
bers mj_ . 

In returning to the notations of equation (21-3)^ the quantities 



1 Ri 


are proportional to constant numbers . 


p 7 ^ 

In other words, 1“^ , and 




are constant. Following this, the functions a(i) and p(t) are 
connected by a unique differential equation. Obviously one may assume 

I = A. Then A^ is a linear function of time and UqA is constant. 

Considering equation (19 -I), which we recall here 




(19-1) 


we obtain of necessity: 



^te 


^ ^ lOvt 


t 


Sa = 


(21-5) 
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The LoitslaxLsky relation is not verified, so that no total simi- 
larity can prevail; a(i|f) and (3(t) are joined by the relation 


a 



( 21 - 6 ) 


which can be subjected to experimental check. 

A complete discussion of the second and third possibility is fore- 
gone, in favor of the case where Z = A is prescribed. Then (21-5) 
becomes 


y^lra' (t) + a"(\lf) 


fa' (7) 

If 


5a(»|/) 


ki-. 

2-A 


P’(flf) 


fp(lf) 


= 0 


(21-7) 


and is represented schematically by 


)Pp(t ) + 0 (2I-7) 

The discussion then deals with three-dimensional rather than four- 
dimensional space. 

Two cases are possible: 

(1) The point a^^(t) describes a straight line (issuing from 0), 
and point |3j_('vir) remains in the plane perpendicular to this line 
(passing through 0). 

(2) The roles of the two points are permutated. 

The first case is not distinct from that which we have studied. 

In the second case, there exist three constants so that 
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'lfa'('lr) _ _ P’(^) + |p(7) 

kvmj^ ni2 


(21-8) 



m2 




= 0 


They axe the conditions proposed by Sedov. They can be checked 
directly by experiment; m2^ cannot be zero, because would be = 

and at the same time a(.]/) = which is not possible; mg and m^ 

cannot be zero simultaneously. 

If m^ ^ "the similarity cannot be total, because one would have 
Uq 2^5 = c"'"® and by (19-I) 

}? = 4vt Uo2t5/2 • C*® , c“ 


which is incompatible with the second equation (21-8). 

Following this, uq and A are solutions of the system of differ- 
ential equations : 


+ “2 

2^ 




= -lOv-^ 




u^ 


= 0 


( 21 - 9 ) 


a(t) verifies the second-order differential equation: 


a 




mg 


4vmi 


'Jf) a' + 5a = 0 


( 21 - 10 ) 
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after which p is determined by the first-order differential equation 


P' + % = 


ITI7 


2vm 


4a* 


(21-11) 


The system (21-9) can be reduced to qimdratiores . If the solution 
cannot be written in finite terms, it still is desirable to integrate 
as far as possible. 

When the first eq\iation (21-9) is set in the form 


R 

& 


1 - k 
10 V 



+ k 


( 21 - 12 ) 


where Rq and k are constants replacing m^^, m 2 , and m^, it is 

found that 


A = 


ctepl^k 


R 



(21-13) 


Whence it is deduced that 


$ = S - 1% L . 

dt R - (1 - k)^l r] 


(21-14) 


The problem is thus reduced to a quadrature. By (21-13), 
follows from R, and Uq = 


The equation (21-10) generalizes equation (20-4) to a certain extent, 
but it is less simple and will not be discussed. 


Other types of solutions satisfying a similarity hypothesis can be 
examined, such as the case where A is very small, that is, where the 
dissipation due to viscosity is low. In that case a length other than 


f (r)dr . 


The discussion 


A must be chosen for length I, such as L = 
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proceeds on the assumption that the viscosity terms of the fundamental 
equation are negligible against the other two. It leaves then 


The calculations are not developed (for further details, consult 
Batchelor's report (ref. 7 )). The sole purpose in this paragraph was to 
give an idea of the various methods that can be applied to solve the 
fundamental equation on the basis of the similarity hypothesis. The 
physical study of the problem is deferred until later. 

22 . Transformation of the fundamental equation in spectral terms: 

The formulas of reciprocity between the correlations and the spec- 
tral functions made it seem interesting to transcribe the fundamental 
equation (17-IO) in spectral terms. (Compare spectral intermedium, 
section 19 •) 


^ = div T' 


(21-15) 


dt 


By ( 11 - 10 ): 


3 


5 





cp(A)dA 


( 22 - 1 ) 


with 



Likewise 


5 


AR = - / e P ^ kScpdA 




(22-2) 
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Hence 


St 


- 2v AR = 







^ + 2vk2cp)dA 


(22-5) 


^ 

This quantity is equal to div T, 

for components ^ocpa* We put: 

a 


— > 

where T 


is the vector having 


div T = 



3 

IT. 

p=i 


The fundamental equation becomes 


— + 2vk^cp = ^ 
St 


(22-4) 


In this equation cp and t are unknown functions of A, t, and 

k^ = ^ Ap^. 

p=l 

The total energy of the fluctuation is 


E 



(it was shown in section 9 that / \|rdA = 0). 


For isotropic turbulences, the spectral function F(k) = 2jtk^cp, 
which does not depend separately on general case 

for cp, but solely on k, is introduced. 
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The factor ij/ can also be expressed by means of the scalar 


T = Uq5(c + 2a) = = ^(cr^)' 


It is actually known that 


div T^= T' + -T = ^(r%)' (22-5) 

r J.2 V ' 


which^ after elementary calculations^ gives 


■Kk) 



sin rk ^p2ijA ' 
rk V J 


dr 


( 22 - 6 ) 


Equation (22-4) reads 


^ + 2vk2F 

8t 


2«k^^|f 


(22-7) 


Putting : 2rtk^ilf 

triple correlations 


= $(k^t)^ where 
and such that 


J is a function connected with 



yields 


— + 2vk^F = ^ (22-8) 

8t 

By the use of the formula (22-5) and the limited development of 
sin rk, it is readily proved that the function $(k) is infinitely 

small as k^, provided that the correlation function c approaches 
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zero faster at infinity than 


1 



But, with the relation 



dk = 0, 


nothing more is known a priori about its form. 


Equation (22-8) is not essentially different from (I7-IO). It 
raises the same difficiilties and is always a single equation between 
two unknowns F and However, it is much easier to make mathemati- 

cally or physically reasonable assumptions about the form of ^ than 
about the form of the triple correlations, as will be shown in the fol- 
lowing two methods, both due to Heisenberg, which take advantage of 
(22-8); the first involves the introduction of the correlations in time 
and will go beyond the pvirpose originally assigned to it, the second 
leads effectively to certain possible forms of the spectral function F. 


25. First theory of Heisenberg (ref. 25): 

The hypotheses to be examined primarily are those which form the 
subject of the last part of Heisenberg's report. We shall utilize the 
ideas and follow the calculations as closely as possible, but change 
the notations and substitute stochastic averages^O for the spatial and 
temporal averages. 

The expressions giving the statistical velocity and pressure of 
homogeneous turbulence read: 




r 


'A 


p=l 


>vpxj3 




P = 





dq(A) (25-1) 


The equations of motion are first written in spectral terms . It 
can be assumed that hj^(A) is derivable on each test with respect to 

8u 

the time, and its derivative, which corresponds to — is designated 

8t 

by b^(>')- 


b*^Bass, reference 5* The calculations are rather difficult. The 
reader who wants to avoid reading them will find a summary of the 
results at the end of paragraph 24. 
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The expressions (23-1) must he introduced in the equations (25-2) 


Su 




— + > ug-i!: + U- 
at p ^axp 


a 


1 ap ^ 

+ V 2_ 


axp ax^ p axo, p axg^ 


a 


(25-2) 


where U is the constant over-all velocity. 

The only term of which the eqiiation merits a more detailed dis- 
cussion is 




Uo — = -1 


/AXA 


A'pe 


■i ^ ^ ^ ^^p^I 


dh^*(A’)dhp(p) 


(25-3) 

Putting pp - A'p = Ap and considering the equation of continuity- 
yields 


= 0 


which had already been used in the form 


n ^ 'll. 

. 17, / e dh,^*(p - A)dhp(p) (23-i^) 

the subscript (p) specifying that the integral is extended over the values 
of p. This quantity is the Fourier transform with respect to A of 


iAp 



dh^*(p - A)dhp(p) 
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By means of some elementary transformations two equations are 
obtained: 


dh^(A) + i 





/ - A)dhp(|i) + iA^U diijA) 

J(b) 




- - vk2dh^(X) 


(23-5) 


^ i I_ V? f - AjdhpC^) 

P ap 


where h^(A) designates the derivative of hj^(A) with respect to A. 

The demonstration of paragraph I7 from equations (25-5) is then 
resumed. The first equation (25-5) is multiplied by dhg;^(A), added up 

with respect to a, i is then changed to -i and the two equations 
obtained are added member by member. 

Putting 


z = -21 ^ A|3 ah^*(A) 
oP 



dh^*(p - A)dhp 


(h) 


(25-6) 


while considering the equation of continuity^ gives 


- y 

St — 


dhjA) 


+ 2vk^ ^ 


dh (A) 
a' ' 


= R(Z) 


(25-7) 


B(Z) denotes the real part of Z. 

After averaging, it is known that when the turbulence is isotropic: 
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Consequently 


A A + 2 vk^A = R(Z) ( 25 - 8 ) 


What was called $ in paragraph 22 is given by the equation 


2 jtk 2 R(Z) = t dA (23-9) 


In order to find a useful expression of Z several hypotheses 
originating with Heisenberg are made. 

First of all, it is known that the statistical functions dh^^CA) 

are of orthogonal increments, that is, that the averages of the pro- 
ducts dh (A)dh *(A') are zero, if the two points A, A' of the wave 

number space are distinct. The first hypothesis is that the statistical 
variables dh^(A), dh^*(A') are not only orthogonal but independent 

(in probability). 

Recalling that the average values dh^(A) are zero and that 
dhoj^^(A) = dh^(-A), it is seen that the average value of 


dh^*(A)dh^*(p - A)dhp(p) 


is zero, except when at the same time 

A = A- p p-A = |ji A=p 

that is, A = |i = 0. If it is ass\jmed that the probabilities of sudden 
jumps of h^(A) are zero, then Z = 0. 

Paragraph 20 dealt with the case where Z = 0. It is certain that, 
in general, Z is not rigorously zero, but it is likely that Z is small 

8F 

compared to the viscosity term, or the term of development In other 

words, the hypothesis of orthogonality is too general, and that of the 
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independence too restrictive, yet it is difficult to formulate an inter- 
mediary hypothesis. But, starting from the independence it is possible 
to obtain, by means of certain interesting transformations, a significant 
expression for Z. This expression is defective since there is contra- 
diction. But it seems likely that it represents a good approximation 
of Z. We shall therefore study it. 


2k. First theory of Heisenberg (continued). Space-time correlations: 

If the average of r(Z) disappears, it is because Z contains an 
odd mimber of factors dh^(7\); the equations of motion are not linear. 

Replace dhp(^i) by the time integral of its derivative dhp([j.) and 

express dhp(p) by means of the equations of motion. Certain terms 

linear in dh are replaced by quadratic terms; products of four fac- 
tors dh appear in Z, and Z is not zero. A kind of technique of 
solving the equations of motion is involved here, which reduces them 
to integral equations (bearing on simple averages). 


If T is a very great positive, number 




"'t+T 


f-T 


dhp(|a,T)dT = / - T)dr 


(24-1) 


If T is sufficiently large it may be conceded that hp(|j:,-T) 

tends toward zero, that is, that the fluid has "started from rest." 
Hence 


dhp(^,t) 


0 


- T)dT 


(24-2) 


To lighten the notations^ h'p((i) is to represent the func- 
tion hp(^i,t - t). Then 

n 


= -2i <3h *(?\) dT 

aB Jo 


dh 




i^*(b - A)dh’p(n) (24-5) 
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and dh'p(^) is eliminated by means of the equations of motion from 
which we removed previously the term of the pressure in dq: 


dh'p(q) = 


-1 




-q)dh’^(q') + 


2 ^ 
h 





q)dh'pt(q') - (vq^ + iqiU)dh'p(q) 


(24-4) 



Nov a nev assumption is made according to which the original assump- 
tion of independence (for t = O), which is certainly verified for the 
great values of remains valid for the whole time interval leaving 
only the averages with four factors 




H^dh^*(A)dh^*(|a - A)dh’ 






2 ^ r r . A)dh'^,*((jL- - n)dh' ,(q’) (24-5) 

apa'p' ^2 Jo J(h,h') 

The dh being independent in distinct points of the frequency space, 
it is necessary to associate the four points which correspond to the 
foiar dh by pairs, not forgetting that dh^*(A) = dhj^(-'X). Considering 

the transformation of the first integral, the only admissible combinations 
are 


( 1 ) 


st 


A=A-q q'-q=q 

hence q = 0. No contribution to the integral. 

(2)nd 

A=q-q' q-A=q' 
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conditions which are reduced to |j' = |i - A, and leave an integral other 
than zero 



^^dh^*(A)dh'p(A) 


dh^*(^x 


^)dh'^(^ 


A) 


(24-6) 




A = |i' ^l-A = n-^I' 

conditions which are reduced to n' = A, and leave behind the integral 


/ ^^dh^^(A)dh’^(A) dh^*(^ - A)dh-p(^ - A) 

J(\^) 


(24-7) 


After some changes in the denomination of the subscripts, the 
first term of Z assumes the form 


-a^r 

In 


n 


aPr Jo 

^ (24-8) 


Putting p - A = A', while considering the incompressibility condi- 
tion, eqiiation (24-8) becomes 


-2 ) A^A'p 
aP7 


dT dh *(A)dh'g(A) dh^*(A’)<aii'y(A') 

0 J(A') 


(24-9) 


A similar calculation which need not be given in detail gives the 
expression of the second term of Z. Finally 


Z = 2 dT / 

Jo J(A-) 


1 2 , ,2 

k - k' 


k2 + k’2 + 2 ^ ApA 


5'' 6 


ap7 


\^A'pdh^*(A)dh'p(A) dh^*(A')dh'^(A') (24-10) 
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Now the components of the correlation tensor in space and time are 
introduced: 


Sap(l^'t,T) = u^(x,t)up(x + i,t - t) (2A-II) 

The u„ are expressed in spectral terms. Putting 

Ou 


dh *(A)dh’3(X) = dh *(7v,t)dh3(A,t - t) = 0a,R(>^,t,T )dX 

(24-12) 


it is seen that 






Xpip 




(24-15) 


a formula which ties the correlations to the generalized spectral 

tensor Now, Z is expressed by means of 


Z = 2 dA 


Hoo 

dT 

0 u 


k*=^ - k'‘ 


(^')k2 + k'2 + 2 YZ ^6^' 


y~ Vp9^p(A,t,T)G„^(A',t,T)dA’ 

ap7 


(24-14) 


For easier reading of this formula, the term under the sign 
the scalar product of the two vectors having for components 


represents 


^ ^ ^ ^ ^ ) and ^ 9 ^)^^ 


The real part of Z is obtained by transforming the 0^j3 
correlations S^p by (24-13) and taking the real parts of the 
exponentials . 


into 
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All these calculations do not necessarily assume isotropic turbulence 
But in this particular case they are a little simplified and will be term! 
nated. If there is isotropy and incompressibility: 


= A(k,t,r) 




k2 i 


(24-15) 


A 


1 V 


■"oca 


being 


A being real^ 


a real quantity to be interpreted right away. 


Z is real and equal to R(z) . Therefore 


pco n k2 - k'2 




5^ 6 


V 

z_ 


apr 


A^A'pA(k,t,T)A(k' ,t,r) 




A'g^'y 

V 



dA' 


(24-16) 

To teminate the calculation^ the last factors are developed^ and 
with the vector A taken as support of the third axis one passes to 
polar coordinates. The angle 0 of vectors A and A' is introduced^ 
defined by 



5 = kk* cos 0 


Finally, one puts cos 0 = ^ and introduces the elementary function 


J(k,k') 



- k'2)kJk’5 r at 

U~1 k2 + k'2 + 2kk'^ 


= rt ^k 


- k 


.2 


k^ + k' 


- 2k' 

3 


2^.2 


n 

2 


(k2 + k'^')(k2 

kk' 



k 4- k' 
k - k' 


(24-17) 
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Suirimary and conclusion .- By the successive assumptions of the pre- 
sent paragraph, the fundamental equation (22-8) reads 


The function j(k,k') given by (24-17) changes sign when k and 
k* are permuted. The integral of the second member of (24-18) with 
respect to k, from 0 to 00 ^ is therefore zero and (24-18) contains, 
as it should be, the known equation (l9-2): 


8 t 

as it follows from formulas (15-4) and ( 15 -IO). 

The function A(k,t,T) is tied to the space-time correlations: 



dT I j(k,k* )A(k,t,r)A(k* ,t,r)dk* (24-l8) 


loo 


= u^(x,t)up(x + |,t - t) 


by the relation 





P A(k,t,T)dA 


If the turbiilence is isotropic it gives, more simply: 




I' 


1 00 


k2 sin rk^ (ic,t,T)dJc (24-19) 


rk 


and conversely 




being naturally a function of r alone- 
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If T = 0, the general tensor So^p(^,t,T) is reduced to the corre-^ 
lation tensor in space Rj^p(|,t) and ^ becomes E ^oa ~ 


becomes 


F(k,t) = l4-jtk2A(k,t,0) 


(2l^-2l) 


is verified whereby the formulas (24-19) and (24-20) contain the formu- 
las (l4-4) and (l4-5). Therefore, (24-18) is a nonlinear integral equa- 
tion, in which appear at the same time A(k,t,T) and A(k,t,0) = F(k,t) 
and which makes it possible in principle to determine the general 
correlations, in time and space, of isotropic turbulence by means of 
their spectral function A(k,t,T). 


25 . Second theory of Heisenberg: 

This "second theory" came first in Heisenberg’s report (ref. 25)- 
By means of postulates of physical character, it effectively enables the 
fom of the spectral function F(k) to be determined in certain conditions. 
But the formula set up still contains an indeterminate constant. It seems 
that Heisenberg’s idea was to caxry the expression obtained for F in the 
nonlinear equations deduced from the equations of motion, which would have 
made it possible to calculate the constant. The theory of sections 25 
and 24 does not exactly come up to these expectations . It results in a 
nonlinear integral equation in which not the spectral function F(k,t) but 
its space-time extension A(k,t,r) is involved. If A(k,t,0) = F(k,t) 
are known, by the theory which is to be discussed, it seems theoretically 
possible to deduce from it A(k,t,r) for every time interval t. But 
the mathematical problem thus posed appears difficult and does not fit 
into any class of known problems, and no attempt is made to solve it. 

The dissipation of energy per unit mass in the ensemble of the spec- 
trum is, as is known 



€ = 2v 


(25-1) 
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The dissipated energy per unit mass in the spectral range 
> A < k , that is , by the eddies of "dimensions greater than i 

and of relatively large scale are designated by If k — >00, 

tends toward e, and if k >0, ej^. tends toward zero. Heisenberg 
reduced to a form similar to e by introducing an apparent coef- 
ficient of viscosity so that 

+ Vk) J k'2F(k')dk' (25-2) 

Vk must tend toward zero when k — ><». 

Now this first hypothesis must be perfected in such a way that Vk 
assumes a usable form. 


Since lim Vk = 0, Vk can be represented by an integral 
k — »oo „ 

Poo 


What about the function under the sign 


? It is assumed that this 


function depends on the values k' ranging between k and 00 directly 
and by means of F(k'), but that it does not depend on k. In that case, 
it can be determined by dimensional consideration. Its product by dk' 
having the dimensions of a viscosity v, it is necessarily of the form 


C ISEn. where 
1 h'5 


C is a purely numerical constant. 


Hence 



k'2F(k' )dk' 


(25-5) 


ek is the sum of two terms. The first corresponds to the part of 
the energy of the fluctuations of the spectral range (0,k) which is 
transformed directly into heat. The second is the part of that energy 
which, before being transformed into heat, serves, first, to maintain 
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the energy of the fluctuations of the spectral range (k,oo) of the "small 
eddies . " 

Now, the fundamental equation (22-8) can he decomposed precisely in 
this fashion. Integration member by member from 0 to k yields 


8t 


F(k')dk' + 2v 


k'%(k' )dk' 


5(k’ )dk' 


(25-4) 


The first member manifests, first, the energy loss of the fluctua- 
tions of the range (0,k), then the portion of this energy dissipated 
directly by viscosity. The second member is therefore the portion which 
is transformed into energy of the fluctuations of the range (k,oo), and 
the formula (25-3) gives an expression of the function f. Ultimately 


8t 


F(k' )dk' 


+ 2 ( V + C 



Pk 

I k'%(k')dk' = 0 

Jo 


(25-5) 


This is Heisenberg's fundamental equation^^. For the present, it 
is not discussed in its most general form, but only in a special case. 

Heisenberg *s hypotheses are incompatible with the concept of abso- 
lute steadiness, because (25-5) has no solution independent of t. In 
chapter V, it will be shown that Heisenberg’s equation can be applied 
to the problem of spontaneous decay of turbulence behind a grid, an 
essentially unsteady problem. But, an assumption of partial steadiness 
can be made as will be shown in chapter IV. For the values of k 
greater than a fixed approximate value Rq, it can be assumed that the 

energy dissipated by the large eddies (wave numbers lower than k) is a 
constant. Only the manner of distribution between the turbulent fluc- 
tuations of the small eddies (wave numbers higher than k) and the 


^^It appears in this form in the second of Heisenberg’s reports 
(reference 24 ), but the particular case in question is treated in his 
first report (reference 23 )• 
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dissipation by viscosity (transformation into heat, molecular fluctua- 
tions) depend on k, the proportion dissipated as heat increasing with 
k according to formula ( 25 - 5 ) • 

It will be shown (in chapter IV) that this assumption of partial 
steadiness is equivalent to an assumption of statistical equilibrium or 
similarity, which is the more accurate the higher the Reynolds nijmber of 
the turbulence. On these premises the hypothesis is discussed from the 
mathematical point of view. Equation (25-5) must be replaced by 


2 V + C 


dkjj k’2F(k’)dk' = C'te = £o if k > (25-6) 


To solve it, simply put 


9(k) = I k'2F(k')dk' F(k) = ' (k) (25-7) 

lo k2 


The function 0(k) cancels out for k = 0 and verifies the equation 


V + 



f e’(k') 

V.5 


dk' 


€Q 

20 (k) 


(25-8) 


which by differentiation becomes the elementary differential equation 


LL = 

0^ \^0/ k5 


(25-9) 


Putting h 5 
is seen that 


4 fo 

5 20 


and introducing an integration constant 


ki 


it 


0(k) 





4 

5 


(25-10) 
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and finally 



(25-11) 


This formula defines the spectrum of isotropic turbulence for k > 
in terms of an indeterminate parameter and by means of Heisenberg *s 

assumptions^ themselves suggested by Weizsacker’s hypothesis of similarity 
(reference ^2), which will be examined later. 

Discussion . - If k is small against kp, F(k) is reasonably pro- 

portional to k If k is large ^ F(k) is proportional to k""^ (fig- 
ure 7)« These two laws have been indicated by Heisenberg. The first is 
of special interest and used for deriving the corresponding form of the 
correlation functions . 

The formula (25-11) has been established only when k > kQ and^ 

actually^ has no physical significance if k — >0, because F(k) cannot 
become infinite. It is supposed here, in first approximation, that, if 
k < kQ, the true form of F(k) is F(k) = 0, and, moreover, that k 

decidedly is smaller than kp. To interpret this hypothesis mathematically 
is the same as replacing kp by oo in ( 25 -H); hence, one uses the 
physical limiting conditions in which kp is infinite (very high Rey- 
nolds numbers ) . 

Putting 


F(k) = 0 


if 


k < 



( 25 - 12 ) 


if k > kQ 

/ 


the correlation function R(r) can be computed. 
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According to (l4-4): 


5 Ooo _ 5^ 

R(r) = 2F(ko)kn5 | § l n . 3dk (25-13) 

* 


This integral is convergent when kQ > 0; it is not derivable twice 

■•ander the integration sign. An attempt is made to find the form of func- 
tion R(r) for small values of r. It is shown that there exists a num- 
ber a such that 


R(r) 




A being a constant. 

When r — >0, it is easily verified that R(r) has a finite limit; 
equal to 

^ Ooo _ ^ 

3uq 2 = 2F(ko)kQ5 / k 3dk = 3F(ko)ko (25-14) 

J^O . 


Forming 

3uq 2 - R(r) = 2F(ko)ko5 ^1 - ^dk (25-15) 


and making the change of the variable rk = x in the integral, leaves 



(25-16) 
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Posting 


I = 1 - 


Sin x\ 5 


5 

X " dx 


The integral 




differs from I hy infinitesimally small 


quantities with r. Hence, in first approximation when r is not too 
great, 


R(r) = - Ar5| 


(25-17) 


with 


A = = 0.801 


2 

The exponent looked for, a, is equal to — . This formiila is impor- 

5 

tant, and we shall give later on demonstrations for it, of wholly differ- 
ent appearances. The formula that gives the correlation function f(r) 

2 

is similar to (25-17) • H again is a "law in r^. " 

It is a limiting law for infinite kj^. If k^^ is finite, form- 
ula (25 -11) yields for R(r) a law of the ordinary type 


R(r) = 5 uq^ 



which proceeds^ for r = 0, with a horizontal tangent and a finite radius 
of curvature equal to A- When this radius of curvature tends 

toward zero and the limiting curve has the form (25-17) (figure 8), 
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Note.- Comparison of Heisenberg's equations (25-5) (24-18). 

Equation (25-5) is an equation in F, which, theoretically, identifies 
F and can be used to reduce it to a (complicated) differential equa- 
tion. The calculation was made in the "steady" state. It leaves, in the 
expression of F, the constant C, which equation (25-5) can obviously 
not determine. If (24-18) had been an equation in F too, it could have 
been used to determine C . But the procedure that results in this equa- 
tion and whose original pixrpose is to eliminate the triple correlations, 
meets the expectations only partially, because it replaces in some way 
the triple correlations by the correlations of space and of time . Eq\ia- 
tion (24-18) serves therefore, as far as the mathematical problem is 
accessible, to compute the form of the space-time relations, once the 
space correlations are determined by (25-5)* But it does not appear that 
the numerical solution of (24-18) has ever been attempted up to now. 


CHAPTER IV 

THEORY OF LOCAL ISOTROPY AND STATISTICAL EQUILIBRIUM 


26 . Introduction : 

The importance attached to the study of isotropic turbulence was, up 
to the last few years, justified by considerations of, seemingly, essen- 
tially practical value. 

The development of turbulence toward isotropy had in its favor a 
rather feeble theoretical argument, namely, that, far from the sources of 
turbulence, the Reynolds stresses could be assumed zero (see section 6), 
which is a condition necessary, but far from being sufficient, for 
isotropy. 

Batchelor (reference 12) proved in 1948 that this argument is of no 
value; the elements of the anisotropy of turbulence during its formation 
could be recovered in the final phase of its evolution. 

But, from the experimental point of view, the isotropy of turbulence 
behind a grid seems well established, by numerous experimental works, as 
-vfili he shown later. But none of these evolved a general concept; it 
seemed that a particular case was studied here solely by reason of its 
simplicity . 

The theories to be examined here form, by way of contrast, an vinusual 
ensemble, not only because of the successful experimental verifications 
but particularly because of the scope in which they contribute to further 
studies of turbulence. They bring to light the characters of turbulent 
motion which are of general significance. 
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The most outstanding resiiLt is perhaps the better understanding of 
the role played by the viscosity in turbulence . The viscosity is 
involved at the birth of turb\aence and at its death; at its birth it 
is the condition necessary for the creation of the rotation and the 
diffusion of large turbulent masses in the fluid; it intervenes at its 
death (decay) since the turbulent energy is ultimately transformed 
into heat; but the complexity of the turb\ilent motion stems precisely 
from the fact that the viscosity plays a negligible part during the 
major part of the transformations of turbulent energy. 

The fundamental outline is as follows : relatively large turbulent 

masses with relatively low velocity gradients are formed in the creative 
regions of turbulence. The sense of the word relative can be defined by 

introducing the Reynolds number — — , where is the mean kinetic 

V 2 

energy and Lq the mean quantity of these t\rrbulent masses.^ In order 
to have turbiilence in the proper sense of the word^ the ratio 


PUp ^ 

pv^ 

Lq 


V 


must be high. If it is small, we find oiorselves at once in what Batchelor 
has called the final period of txirbulence, in which the rotation decays, 
so to speak, on the spot. The forces of inertia play a secondary part. 

The flow is quas i -laminar . This is the phase through which all turbulence 
passes . It has been experimentally and theoretically studied by Batchelor 
and Townsend (references 10 , 12 ), and we shall come back to it later on. 

If the Reynolds number is high, the vortices undergo a double evolu- 
tion, according to Heisenberg's representation (reference 23); they get 
bigger on the one hand and smaller on the other, that is, their spatial 
dimensions increase on the one hand, and break up into smaller vortices 
on the other. Leaving aside the first aspect of this evolution, for the 
time being, the theory of similarity studies the development of the energy 
passed on to the small vortices; the latter undergo a development similar 
to that of the generating vortices, they pass on their energy to smaller 
vortices and so forth up to the moment where the size of the vortices is 
such that, their Reynolds n\amber is small. A bigger and bigger part of the 
energy is then changed into heat. This somewhat vague representation is 
defined a little more accurately fiorther on, in particular, by the analysis 
of the rotational development. 

■^^This rather vague expression is specified to some extent by the 
study of dynamic equilibrium (chapter V). 
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It is reasonable to think - and this assumption will be backed up 
by its consequences - that the development of vortices is the more rapid 
the smaller their size. On the other hand, the statistical character of 
the velocity of the fluctuations, which is prescribed by experience, is 
assumed. 

Hence it may be imagined that, when the energy has been transmitted 
to vortices sufficiently smaller than the initial vortices, the structure 
of these vortices is., on the average, independent of that of the generating 
vortices , and that the relatively slow development of the latter has no 
effect on the development of the smaller size vortices. The latter, there- 
fore, find themselves in a quasi -balanced state, one determining factor of 
which is the energy per unit mass and time which is supplied by the bigger 
to the smaller vortices. 


27. Definition of local homogeneity and local isotropy: 


At first we shall follow Kolmogorof f ' s method (reference 31. For 
more details see Batchelor, reference 6) which seems the most direct as 
well as the easiest, at least, to begin with. Kolmogorof f introduces 
quantities which can be completely defined in a limited spatial domain G 
and a limited temporal field T. Assume that M and M' are two points 
of G, that t and t' are two instants of t; we call r the dis- 
tance MM ' , and put r = t - t ' . 


The difference w of the velocities at M and M' 
and t' is introduced 


at instants t 


Wi = \ij^(M,t) - Ujj^(M',t') 


(27-1) 


and defined by its law of probability in the space-time domain (G,T). 

One notices immediately that w is little dependent on the spectral 
components of the velocity which develop slowly in the spatial domain G 
and in the time interval T. 

By its definition w depends a priori on M and t, as well as on 
M ' and t ' . 

Kolmogoroff defines what he calls the local homogeneity and local 
isotropy. A condensed version of these definitions reads as follows: 

The turbulence is said to be locally homogeneous when the laws of 
probability for the statistical function w depend, in the domain (G,T), 
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only on t and the vector MM' (actually^ instead of this vector, the 
vector ^ defined later on is involved) . 

The turbulence is termed locally isotropic when it is locally homo- 
geneous and when these probability laws depend neither on the orientation 
nor on the sense of the coordinate axes. 

According to the foregoing, it may be assumed that in every turbu- 
lence having a rather high Reynolds number, there exists a spatial 
domain G which is small compared to the size of the big vortices, and 
a time interval T, short compared to the time necessary for their energy 

to change appreciably, and such that the turbulence is locally homogeneous 
and isotropic in the domain (G,T). 

But these definitions can have no physical significance unless the 
reference system used by the observer is itself physically coordinated 
with the turbulence. The domain G must, in some fashion, follow the 
turbulence with the mean flow velocity. Kolmogoroff, in his original 
report (ref. 3I) defines the velocity w about a reference point (MQ^tQ) 

in terms of the components of the vector MqM and of t = t - tQ 

by the formulas 


= u^(M,t) - Ui^Mo,tQ^ (27-2) 

then he introduces a vector e having as component 


e^ = - TUi^MQ,to) (27-3) 


This is a random vector whose physical significance is evident in 
figure 9. Since the essential point of the vector Tu|MQ,tQ^ arises from 

the "mean velocity" in the domain G, the use of the coordinates e^^ 

is tantamount to following the turbulence in its entire motion. 

The correlation measurements in time relative to the laws of proba- 
bility of w are few and far between. Some measurements have been made 
in water and in air, but none in the wind tunnel, as far as we know. 
Obviously, in such measurements, time-space intervals must be combined 
and it must be ascertained whether correlations close to 1 can be 

obtained for a finite value of the ratio which, if it exists, define. 

T 

the mean velocity at the scale involved. 
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But, at the present state of research the study is limited to simul- 
taneous quantities for defining the statistical values w. As in ordinary, 
isotropic turbulence, it is usually limited to the moments of the second 
and third order. The tensor of the moments of the second order 


Wi^M,Mo)wj(M’,Mo) 


simxilteineously introduces three points, but can be expressed, on account 
of the isotropy, homogeneity and incompressibility, as a sum of moments 

relative to the couples ^M,MqJ, 




Wi (M ’ ,Mq) w j (m ' ,Mq) - Wi (^M,M 'j w j ,M ') 


(27-^) 


In this case, the tensor 


B 


ij 





plays a role parallel to that of the tensor Rj^j in isotropic t\irbulence, 
and by the same argument it is apparent that the isotropy leads to writing 






(27-5) 


where 


Sdd^'') = (""'d - ""d)^ 


B (r) =■ (n' - u ^2 

nn' ' V ^ W 
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and representing the velocity components along MM' or normal 

to MM' . 


A total moment can be defined too: 


B(r) = ^ ^ (u'l - Ui)2 

The incompressibility involves the relation 

Bnn = Bdd + I (27-6) 


equivalent to Karman^s relation in isotropic turbulence (9-7) • Other 
relations can be established^ in particular, the expression of the 
dissipation of energy by viscosity 


e 



. L: 


(27-7) 


Likewise, the moments of the third order can be expressed by the 
single fiuiction 


Bddd 




3 


28. Similarity hypotheses. Statistical equilibrium; 

The definition of these local quantities enables the notion of sta- 
tistical equilibrium of small vortices to be defined. According to our 
hypotheses, this equilibrium is not dependent on the particular charac- 
teristics of the flow. The only quantities which can be of influence 
are : 

(1) The energy supplied per unit mass and time to the small vortices 

by the larger ones; 

(2) The factor governing the kinematic viscosity. 
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If the fluid is assumed incompressible, the specific mass can be 
discounted. 

Hence Kolmogorof f ' s postulate of similarity: in a domain where the 

turbulence is locally homogeneous and isotropic, the laws of probability 
depend only on r, e, and v. 

This postiilate justifies the dimensional analysis which results in 
setting 




1 

In this formula, is a universal function, and I = 

represents a length, called the local scale of turbulence by Kolmogoroff . 
This length is much smaller than the length of dissipation A of the 
isotropic turbulence. It is, in effect (compare (l4-ll)) 




whence 




u^A 

where the Reynolds number R-v = must be at least some tens if the 

—A ^ 

conditions of similarity are to be realized. 

The nondimens ional statistical quantities must be universal func- 
tions of — . Consequently 
I 
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(coefficient of dissymmetry) 




(coefficient of flatness) 




'^(t) 



(28-2) 


and similar quantities Aj^, with the normal velocity component 


They are \miversal fimctions in a certain domain Z « r « L, in 
which L is a length characterizing the size of G. It is to he noted 
that, accorulng to the hypotheses, I and L can he slowly variable 
functions of time. The precise definition of the range of validity of 
the formulas is ticklish, hut in any case, if r is made =0, the 
values of 


S^( 0 ) 






2 


and 


*d(0) 



(28-3) 


must be universal constants^ which affords a first confirmation of the 
similarity hypothesis. 

Townsend (ref. 57 ^ 58 )^ sit Cambridge^ developed a method of measiaring 
S( 0 ) and A( 0 ). The measurements were first made in isotropic turbulence 
produced by a grid. These coefficients are constant over a wide range of 
Reynolds mmibers (ratios varying from 1 to 5) sit varying distances 
from the grid (likewise in a ratio near to 5)^ that is, at different 
moments of the development of the turbulence. 

Townsend quotes the value 5-^9 - 0 . 0 ^ for A and -O.38 for S, with 
a little lower accuracy. It should be recalled that for Gaussian distri- 
butions the values of A and S are, respectively, 3 0 * 
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To check the theory of locally isotropic turbiilence, by which the 
structure of the small vortices is independent of the nature of the large 
vortices and the formation of turbulence , Townsend made also measurements 
in the wake of a cylinder (refs. 39, 40, and 4l), that is, in turbulence 
neither homogeneous nor isotropic (in the original sense). He determined 
thus S(3^(0), 3^(0), A^(0) and An(0) over the entire extent of the 

wake, at distances of more than 80 diameters from the cylinder. The 
values obtained are the same, with sufficient accuracy, as those obtained 
for isotropic turbulence behind a grid. Furthermore, the measurement of 
the quantities 



Uj_ being parallel to the overall velocity (axis of xq), shows that, as 

for isotropic turbulence and throughout the entire wake, these quantities 
satisfy the equations of isotropy 



By the same argument, the expression 


15 

±2v ^ , as for 

2 \bxj 

(27-7) may lead to presume. 


isotropic turbulence, and as the equations (27- 


is 

and 


These formulas are very important because they show the possibility 
of defining the dissipation of energy in any tirrbulent flow by measure- 
ments bearing on a single velocity component. 


29. Case of high Reynolds numbers (reference 6); 

The study of the statistical equilibri\im of vortices can be extended 
further when the Reynolds number is high enough so that the viscous dis- 
sipation in the largest vortices of the statisticaJ. equilibrium is rela- 
tively low; it can be assumed then that the statistical equilibriimi of 
these large vortices is not dependent on v. 
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This is Koimogorof f ' s second hypothesis: 

If I « r « the laws of probability regarding w are then 
dependent on e and r only. 

This implies that a domain G can be fo\md whose size L is large 
compared to I, that is, that the Reynolds number is very high. 

Dimensional analysis permits then the definition of the form of the 
moments of the statistical variable w. In fact it is readily apparent 


C being a universal constant (figure 10 ) . 

The moment must have a similar form. The equation of continuity 


that, for dependent on v 



must be of the form 



whence 


2 


B^^(r) = C(€r)5 


(29-1) 



On the other hand, when 


r 


is smeill 


2 


®nn(^) = |®"nn(0)^^ = 


(29-2) 


Likewise, one must have 


BdddCr) = X cr 
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But, if r is small 


Bddd(^) = X r5 


(29-5) 


^3 

2 


because in that case, ^ddd^^^ f ®dd(^) 

Since I depends on v, the nondimensional quantities such as S^f— 


in first approximation. 

/r\ 


I' 

must be constant, for values of r much greater than Z. On the other 
hand, Sjj^(O) is also a universal constant, and the experiment seems to 

indicate that these two constants are identical. 


To tell the truth, there is only one series of measurements of B^(r) 

made by Townsend in a region where the turbulence was isotropic, by 
measuring the triple correlation function 


c r) = 

3 

^0 


The relations 


^u - u'p = -6uo^c(r) (u - u’)^ = Uq^ [I - f(r)j (29-^) 


are utilized. 

The experimental points are placed near = -O.38 for values of 

— of the order of kOl , The absolute value of S decreases slowly 
I 

afterward. 

Instead of arguing about the correlation functions, one may just as 
well apply the dimensional analysis to the spectral fimctions F(k) for 
values of k which are sufficiently large without being too large, how- 
ever. Since F(k) must be dependent only on k and e, it is found that 

2 _5 

F(k) must be proportional to e^k 5 in the domain under consideration. 
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This is the spectral law in k 5 spoken of earlier during the discussion 
of Heisenberg *s theory (section 25)- 

Using the detailed hypotheses of Heisenberg, we have then shown that 

2 

there is equivalence between the law of correlation in r^ and the spec- 

.5 

tral law in k 5. It is now seen that this equivalence was inevitable, 
and that it is not the result of the particular form of the law of energy 
transfer in the spectrum, but solely due to the existence of a general 
similarity hypothesis among the hypotheses the utilization of which was 
unimportant . 

By using the equations of motion the value of can be defined 

and a relation established between the constants C and S, which is sus- 
ceptible to experimental verifications . 

The equation corresponding to that of Karm^ (l8-5) reads here 


-4e = B* 


ddd 


+ 







(29-5) 


The time factor does not appear, because the hypothesis of local 

d 2 

homogeneity contains that of steadiness with time, so that — uq must 

dt 

he taken equal to - ^ . 

3 

The equation (29-5) multiplied by is integrated in the form 


6vB' 


dd 


B 


ddd 


(29-6) 


If r is small, can be neglected, so that one finds again 

the relation 
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For the values of r in the interval I « r « 

over because is proportional to r and 

Hence 

®ddd = - 

5 

But in (29-3) it was seen that in first approximation: 


L, B 


ddd 


prevails 


B'dd ^ 


1 

5 


%dd.(^) = Sd(0) |B(3^(i(r)j2 _ 


1 

Sa(0)c2er 


Conseq^uently 


C 



(29-7) 


Therefore, C can be calculated, if the value of S^(0) is known 
and, particularly, if the experimental value S^(0) = -O .38 is assumed. 
In that case 


Bdd(^) 




2 

5 


Conversely, the experimental study of the function B^^(r) enables 

C to be meas\ired and the obtained value to be compared with the value 
computed from S^(0). 


50 . Validity of the similarity laws: 

All second-order moment measurements have been made in isotropic 
turbulence. Hence it is usefiiL to examine what becomes of the theory of 
local isotropy if applied to isotropic turbulence proper. 
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The first point is the determination of the space -time range of 
validity. The discussion, necessarily rather vague, has been made by 
Batchelor (ref. 6). It is briefly summarized here: 


From the spatial point of view, isotropy appears to be assured for 


dimensions of the order of magnitude of the correlation length 





The difficulty arises from the estimation of the time interval. An 
evaluation of the characteristic period of the vortices indicates that 
the ratio between the period of the largest and that of the smallest 

vortices is of the order of — , where A is the length of dissipation. 

A 

This ratio must, therefore, be great. Experience indicates an order of 

magnitude of in what Batchelor calls the "initial phase of decay of 

A 

turbiaence." (It is to be defined in chapter V.) For the present, it 
is simply stated that it is produced by the turbulence behind a grid of 

mesh size M, when the Reynolds number ^ of "the grid (U = mean 


velocity) is sufficiently high (1,000 to 300,000) and the meas\irements 
are made at a distance from the grid of the order of 1,000 M. One finds 


T £> 

then that — is of the order of — , where ^ = 

A 10 V 

Reynolds number of turbulence. During the phase in question, R^ is 


UqA 


may be called the 


constant and eqiial to 
constitution of the grid. 


— , where 

A 


A is a number dependent on the 


For values of R, not quite up to 200, which still gives only values 
—A 

of ratio — of the order of 20, Townsend has measured the correlation 
A 

function g(r). Since 


1 - g(r) . 

2Uq2 


1 - g(r) 
2 

r3 


must be constant. 
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The curves plotted for various values of ^ and ^ have similar 

shapes. They present a rather flat maximum between r = O.IL and 
r = 0.5L^ from which the value of the constant C (section 29) can be 
deduced. The average of the obtained values is 1.55> with discrepancies 
of less than 3 percent. Considering the accuracy of the measurements of 
Sd( 0 ), the agreement with the value 1.64 deduced from S(j^( 0 ) is 

sufficient . 

Other measurements, made at a lower Reynolds number, yield C = 1.33^ 
while those by Dryden produced C = I.50. In conclusion, although the 
range of a validity of the law is too restricted for the aspect of the 
curves to be convincing, the agreement of the values of C is quite 
satisfactory. 


The lower limit of validity of the law of local isotropy in these 
experiments is of the order of magnitude of X, that is, well superior to 

I, since — = is here of the order of 20 to 25. 

z 


It seems that the Reynolds numbers used are still too low for the 
field of application of Kolmogorof f ' s second hypothesis to be extended 
very far. The first hypothesis, less restrictive, is much better verified. 
Townsend's correlation measurements show, in fact, that the curves repre- 
senting g(r) are when r « X functions of ^ only, at different dis- 

A 

tances from the grid. The similarity is extended to values of ^ higher 

than unity for the highest values of where the second hypothesis 

possesses a domain of validity. This validity is tied to the fact that 
is constant at the beginning of the turbulent development. (Compare 

chapter V.) This information on the development of turbulence can, 
obviously, not be deduced from the theory of local isotropy which ass\mes 
e constant. 


On the other hand, the absence of similarity in the correlation 

curves at high values of - suggests that the range of application of 

A 

the similarity law is very restricted, and also belief in the fact that 
the ratio of correlation length L to dissipation length A varies in 
reasonable fashion during the evolution of turbulence, even when R^ 

L B.A 

remains reasonably constant, because the relation 7 = — is only a 

X 10 

first approximation. 


In reality, the domain of statistical equilibriimi is extremely 
extended, but this does not appear in the structure of the spatial corre- 
lation tensor. For a clear definition of the individ\xal behavior of each 


NACA TM 1377 


121 


scale of vortex dimensions the introduction of the spectral tensor is 
necessary. Heisenberg *s theory is therefore found capable of exceeding 
the stage reached by the theory of local isotropy, and it likewise 
permits interpretation of the different stages of development of turbu- 
lence discussed in paragraph 51^ an^ ia chapter V. 

51. Interpretation of the laws of statistical equilibrium in spectral 
terms - Weizsacker*s and Heisenberg *s theories (refs, k-2 and 25^ 
respectively) : 

Kolmogoroff *s theory assumes a locally steady equilibrium and its 
spatial domain of application is necessarily restricted. However, exami- 
nation of the correlation curves and of the laws of variation of the 
correlation length L during the decay of isotropic turbulence behind a 
grid leads one to believe that the scope of similarity will be extended 
farther than anticipated. 

But this notion of similarity can not appear clearly unless the 
influence of the different vortices is separated. This means changing 
from the correlation functions to spectral functions. 

It will be remembered (compare sectio ns 11 and 12) that to the 
spatial correlation tensor R(xp(^) = there corresponds the spectral 

tensor T(xp(^) defined by 


If there is isotropy, the spectral tensor depends on a single spectral 
function F(k), just as the spatial correlation tensor depends on a single 
function R(r) for example. 

To the properties of H(r) for small values of r correspond those 
of F(k) for great values of k, as results from formula (l4-4) which 
links F to R, and which is recalled here: 


5 




R(r) = 2 


rk 
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It has been shown (section 29) that to the law of correlation in 

2 

r^ for R(r), applicable in the domain I « r « L, and res\ilting from 
Kolraogorof f ' s second hypothesis, there corresponds a spectral law of the 
form 

2 .1 

F(k) = Ae^k 5 


where A is a universal constant. This law is valid in the domain 
kQ< k< kg and the numbers Zkp, LkQ are of the order of unity. 

In the most general case of Kolmogorof f ' s first hypothesis, one can 
only put : ^ 

F(k) = Fo«(i) 

for k > kQ, Fq and k^ being fionctions of e and v only and $ a 
universal function. 

Pk 

Putting J ^(k')dk’ = T(t) in (22-8), gives 


Bt 


F(k')dk' = -T(k) 


V 


Pk 

2k'2F(k')dk’ 

Jo 


(31-1) 


The function T(k) represents the transfer of energy by turbulence 
of the large eddies to smaller eddies. 


In the case in which the dissipation by viscosity is negligible, 
the conditions of statistical equilibrium are written simply T(k) = c 


te , 


5 

5 

A dimensional analysis affords the spectral law in k . In fact, it can 

Bun 

be said that T(k) has the dimensions of averages so that is 
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v^5k^ where is a certain characteristic velocity of vortices having 

k for wave number. Hence one must have V]^5k = c^^, in the domain 

_1 

kg « k« kQ, and v^^ is proportional to k 5. 


For velocity vj^, the quantity 



which, according 


2 

to Kolmogorof f * s theory, is proportional to r^ in the chosen conditions, 
can be assumed. Weizsacker's law is equivalent to Kolmogorof f ^ s law 


because 


F(k), having the dimensions of 



is proportional to 



(Compare section 25-) 


The variation with k of other quantities associated with different 
vortices can also be determined. For example, the characteristic period 

of the vortices is proportional to — a quantity which has the dimen- 

kVj^ 

2 

“ 5 

sions of a time interval, that is, to k . It decreases when, k 
increases, which is in agreement with the hj^othesis according to which 
the characteristic periods of development of the small vortices are 
shorter than those of large size vortices. 

It can likewise he verified that the viscous dissipation, which is 

1 

proportional to k%, varies like k^. Hence, it cannot be neglected for 
great values of k. 

If now the viscosity is no longer neglected, the form of the fimc- 
tion ^(k) must be specified. This is what Heisenberg did who, as 
shown in section 25? bas assimied that T(k) coiild be put in the same 
form as the losses of energy by viscosity, the form of the "tiirbulent 
viscosity" coefficient resulting from similarity considerations . 


The comparison with experiment of the spectral law k-^ has not as 
yet been carried out on flows in which the Reynolds mmiber is high enough 
so that the conclusions are well defined. Nevertheless, its range of 

2 

validity seems more extended than that of the law r^ for the correla- 
tions, because the influence of the various vortex sizes is different due 
to the use of the spectral analysis. 
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Townsend attempted to check the validity of the expression of the 
function T(k) given hy Heisenberg, and according to which F s hould b e 

proportional to k""^ at great values of k. The mean value of 
is proportional to k^F(k)dk. If the law in k"'^ were verified, this 

Jo 

integral should be infinite, which signifies that the integral 



Pk, 


k°F(k)dk should increase considerably together with its upper 


limit k 2 » However, the experiment seems to indicate that this integral 
has a limiting value independent of k 2 , and variable during the decay of 

turbulence. Therefore, the law in k"T cannot be true for the highest 
values of k, which limits the range of validity of Heisenberg's theory. 


CHAPTER V 

DECAY OF THE TURBULENCE BEHIND A GRID 


52 . History: 

The theory of local isotropy shows that the laws of tiirbiilent flows 
are the same, regardless of the origin of turbulence, provided that the 
conditions are "local" and the Reynolds numbers high enoxigh. This chapter 
deals with the particular laws of a turbulent flow of great importance, 
that which is produced in the test chamber of wind tunnels, and in which 
turbulence is largely due to a grid, that is, an obstacle of periodic 
structure. Turbulence arises from the mixing of wakes of grid elements, 
and is dissipated progressively as it moves farther away from the grid. 

For an observer carried along with the velocity of the ensemble, the 
turbulence is therefore subjected to a "decay" as function of the time. 

The laws of this decay axe, naturally, compatible with the local laws of 
Kolmogoroff, but today they are known with a high degree of acc\iracy in a 
much more extended domain, and it can be stated that it is in the study 
of this phenomenon that the theory of turbulence, steadily checked by 
experiment, has made the fastest and most constructive advance. 


17B 


NACA TM 1577 


125 


The basis of this theory was established by G. I. Taylor, who, by 
considerations of similarity, had indicated that, if the Reynolds number 


is high, the turbulent energy 
of the time. 


1 2 


must decrease as t 


-1 


as a function 


The measurements made since by numerous researchers have shown that 
this law is close to experimentaJ. results; nevertheless the e^q)erimental 
curves were rather of the type 




n raiiging between 1 and 2. As the precision of the measurement is low 
as soon as the turbulence intensity becomes low, the Interval to which 
the comparisons refer, is restricted. 

The variation of the length of dissipation A in terms of time is 
tied directly to that of the intensity of turbulence by the relation ( 19 - 1 ) 

1 duQ^ _ lOv 
uq^ dt A ^ 

If Uq varies as t ^ one deduces from it: 


A^ 


n 


But, up to the last few years, the measurements of A offered no 
possibility of defining the value of n. The value of A is, in effect, 
deduced from the curvature at the origin of the curve representing the 
correlation function f(r). The measurements are made with two anemom- 
eters, of which the distance r cannot be reduced indefinitely for 
manifold reasons (geometrical and physical). As a resiilt, the curve f(r) 
is not well known for small values of r. 

However, an examination of the correlation curves seems to indicate 
that their shape remains constant in terms of the distance from the grill, 
during the decay of turbulence. 

All these questions have taken a decisive forward step, as a result of 
A. A. Townsend *s new methods which enable direct measurement, by electrical 
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methods, of the momentary or mean values of the quantities u, — , 

and those of the products of two or three quantities^ as well as the 
direct determination of the statistical frequency curve of a quantity in 
terms of the velocity. 

The use of these methods enabled G. K. Batchelor and A. A. Townsend 
to detemine the law of decay of turbulence produced by a grill with a 
good degree of accuracy over a wide range of values of the Reynolds num- 
bers of the grill, defined by 


p _ UM 


(32-1) 


where M is the width of the mesh, U the overall velocity. 
When the turbulence is isotropic 




dt 


(32-2) 


The direct measiirement of 




affords the evaluation of A 


much more quickly and more accurately than by the procedure employed 
previously, and consequently defines A in temns of time, that is, of 
the distance from the turbiilence -producing grid. 


55. Initial and final phase of tiirbulence (Batchelor and Townsend, refer- 
ences 8 and 11): 

Townsend's measiirements made it possible to distinguish in the decay 
of turbiilence an initial phase and a final phase, separated by a less 
well defined intermediary phase. 

The initial phase, clearly proved by these experiments, corresponds 
to grill distances of less than 100 or I50 mesh widths (variable according 
to the Reynolds number). It is characterized by the fact that the decay of 

turbulent energy can be represented by an expression proportional to t"*^ 
where the exponent n differs from imity by less than 10 percent. It is 

2 -1 
said that Uq decreases as t 
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In the final phase, which corresponds to the limiting state of decay 

_5 

2 2 

of turh\ilence, experiments prove that decreases as t 

The initial phase is also characterized by the fact that the Reynolds 
number of turbiilence, defined by 






(53-1) 


remains constant thro\ighout the period of decay to which it corresponds. 

If, in fact, Uq^ = At"*^, it is seen that that 

Thus remains constant if n = 1, and in this case only. 

To define the structure of the two essential phases of turbulence, 
we now introduce the rotation (Batchelor and Townsend, reference 10) of 
which the relations with the correlation functions and the dissipation 
length A are simple, at least when the turbulence is isotropic. The 
general equations relevant to the rotation are simplified by the Karman- 
Howarth equation. 

The rotation of the velocity, defined by the formulas 


dua 

du2 

du]_ 

du^ 



dx^ 


dxj^ 

CD3 = — 
OXj^ 


- — (35-2) 

OXo 


satisfies two well known equations which are obtained by simple combina- 
tions of the equations of motion 


Su. 

St 


1 ^ 1 dp 

+ 2_ + V Au^ 

I dxj 


P dx. 


(33-3) 
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If (35-5) is derived with respect to since one interchanges 

i and j and subtracts the two equations obtained, the result is 



Tr — 8 o)l5 V — 

V + V £<01 

k oxj^ k oXj^ 


(53-4) 


The two members of (55-4) are multiplied by 2coq, or added with 
respect to with due allowance for incompressibility^ the averages 

2 \ 2 

taken and homogeneous turbulence assumed. Putting leaves 


St 


= 2 2_ ^ 

Ik 


r 

ik 




L 

8xb 




(33-5) 


If, in addition, the turbulence is isotropic, then (section 9): 


0^ = 


% 


= “2^ = 


0^ 


5 _ 


= 5 



In this equation^ the first term of the second member represents 
the increment of the vorticity due to the elongation of the vortex 
filaments^ and the second the decrement due to the influence of the 
viscosity. It is verified that this last term really is negative because 


2o:^ 


8xk^ 


8xk^ 



and =— is zero if the turbulance is homogeneous . 



Batchelor and Townsend meas\rred the three terms of (35-5) directly. 
But it is convenient to give equation (55-5) a slightly different and 
easier-to-interpret aspect. The K^rm^-Howarth equation (l8-5) can be 
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used to evalviate the first member. Introducing the limited developments 
of f, up to the terms in r^ (cf. section 9) 


f(r) = 1 + ^f"(0) + |^f(^)(0) + . . . 


/ ^ 5 c'"(0) 

i(r) = ^ + 

6 


and identifying the terms in gives 


+ ^^<^'"(0) (55-6) 


Now, the derivatives of the functions f and c can be expressed 
in terms of the average values bearing on the velocity derivatives. The 
calculation, as done before for the rotation, gives 


„0%"(0) . -(^ 








,^^'^>(0) = (JT 




If therefore the nondimensional quantities 


(33-7) 






s = -■ 


ySxi 


G = 


2 



2 - - 

■ / V p“ 

/ 8 ui\^ 

2 

(— ) 

f)J 


_w _ 


(33-8) 
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are introduced of which the first has already played a prominent part in 
section 28 (coefficient of dissymmetry, for r = 0), it is seen that 



-a5c"'(0) 


G = = A^f^^^(O) (53-9) 

f "^(0) 


A being, as will be remembered, a "length of dissipation" defined by 
f"(0) =-J^. 

— „ 2 

Since v? = eq,\iation (35-6) can be written 

= _I_(a^)24 _ (35-10) 

6t 5\f3 \ V 


where R-^ = -2— is the Reynolds number of turbulence already defined by 
—A V 

equation (53-2). 

So in turbulence language, equation (35-10) is a consequence of the 
K^ffi^-Howarth equation of turbulence, equivalent to (55-5 ) > superior 

for comparisons with experiments. 

The measurements of the various terms of (33-10) during the initial 
period of decay of turbulence confirm, first, that 


R = A^ = lOvt 


(55-11) 


t being counted from any convenient origin. 
Consequently one may write: 
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and (33-10) takes the form 


G = 

- 2-^ 


S + 


7 


(35-13) 


It has been shown (section 28) that the measurements of the coeffi- 
cient of asymmetry S prove that the latter is independent of and 

close to 0.58. The direct measurements of G are in good agreement with 
the above formula, in which S = O.38, which constitutes a check on the 
qualities of the measurements. 

These results make it possible to evaluate the ratio of the two 
terms of the second member of the equations (53-5) (33-10) • The quo- 

tient of the term representing the dissipation through viscosity to that 
which expresses the inertia effects is equal to 


^ = 1 


60 

7SR^ 


(33-1^) 


or, replacing S, by its experimental value 


1 + 


22 


(35-15) 


This ratio is constant throughout the initial period of decay of 
turbulence. If — is considerably less than unity, the contributions 
_ ^ 

to of the viscosity and of the inertia terms are of comparable 

dt 

importance and the velocity resulting from the variation of the vorticity 
is of an order of magnitude lower than each one of these two terms. This 
is the range of high Reynolds numbers R^, where the theory of statistical 

equilibrium (Kolmogoroff , Weizsacker) applies. For large enough values of 

_5 

the wave number k, the spectral law k ^ can be observed there. The 
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measurements show, in fact, that, in proportion as ^ increases, the 

correlation curves approach the limiting shape with tangent vertical to 
the origin as implied by the formula 


f = 1 



The original parabolical region becomes consistently smaller as 
becomes higher (Figure 11; compare also the end of section 25). 

If approaches unity (Reynolds numbers neither too high nor too 

low), the results deteriorate and new partial or total similarity hypoth- 
eses must be made, resulting in solutions of the type of those studied in 
section 21 from the purely mathematical point of view, and which are not 
repeated here . 

If — is great (low Reynolds numbers), the effect of the triple 
-X 

correlations is negligible against the viscous dissipation. The mechanism 
of the decay of turbulence is mainly due to the viscosity, and the solu- 
tions of section 20 can be invoked. That is what happens in particular in 
the final phase of decay which is discussed below. 


54 . Concepts regarding the structiire of the final phase of turbulence 
(Batchelor and Townsend, reference 12; Batchelor, reference 8): 

It is necessary to predict by theory the law of the final phase of 

_5 

p p 

decay of turbulence, according to which Uq decreases as t and, 
conseq,uently 


= 4vt 


(5^-1) 


When the turbulence is sufficiently atten\xated, the fluctuation is 
slight, and the inertia effects negligible against the effects of vis- 
cosity. As a result, the nonlinear terms in the eq^uations of motion can 
be discounted, and one may simply write 

i- = — + V AUj^ 
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In these conditions the equations of motion can he integrated, the 
velocity components computed in terms of time and the initial conditions, 
and the correlations deduced. This is the program set forth in section 16 
and it cannot be accomplished except in very special cases . Thus the 
time correlations can be calculated just as well as the classical space 
correlations . 

Only the case of isotropic turbulence is analyzed. To neglect the 
inertia terms is to neglect the triple correlations in the "fxmdamental 
equations." Thus one reverts to the conditions of section 20, or t\irns to 
the possibility of expressing the spectral equation in the simplified form 


^ + 2vk2p = 0 {^h-2) 

St 


The general solution of this equation reads 
F(k,t) . 


( 5 ^- 5 ) 


being the instant that makes the start of the final phase . 

One may also use equation (20-1)^ equivalent in correlation terms 
to the spectral equation (5^-2), which^ with the triple correlations 
disregarded, reads 


^ = 2v AR 
St 

of which the general solution (20-9) is a little harder to write. 

But these solutions must have an asymptotic character. It is of 
no interest here, except when t - tQ_ is siafficiently great. F(k,t) 

approaches zero when t - t-^ approaches infinity, as is natural, since 
the turbulence decays progressively. So, when t - t^ is great, the 
only regions of the spectrum continuing to exist are those for which k 
is small, k^^t - t]_j remaining practically finite. It has been seen 
(in section l8) that, for small values of k, F(k) was necessarily of 

the form Ck^, C being a constant independent of the time, and tied 
to the Loi'tsiansky invariant. 
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Consequently, for the large values of t, the spectral function 
takes the form 

F(k,t) = Ck^e"^’^^ (3^-^) 

valid for small values of k, that is, for large-size vortices. 

Example 1 in section 20 gave the corresponding form of the correla- 
tion f\mctions R(r,t) and f(r,t). In particular it is recalled that 

r^ 

f(r,t) = e 


and that 



F(k)dk = - 2 - 

52 v/5 




This is the law of decay that had to be explained. The proof rests 
on the hypotheses of section 20 and section l8, particularly on the fact 
that the Loitsiansky invariant is finite and other than zero. The agree- 
ment of theory and experiment confirms them. 

In this calciilation, it was supposed that the turbulence was isotro- 
pic. On the other hand, only the large size vortices continued to exist 
because the energy contained in the region of the spectrum corresponding 
to substantial values of k becomes negligible. As the constant C is 
independent of time, the shape of the spectriom for small values of k is 
independent of the state of decay of turbulence. The large-size vortices 
contain very little energy, but, in the final phase, the energy of the 
small vortices is dissipated and only the large vortices remain visible 
whose energy, as feeble as it is, has become very great proportionally. 

Now, while the turbtilence is always locally isotropic and exhibits 
a tendency to isotropy for the small vortices (large values of k), it is 
no longer the same for the large -size vortices, in which the geometrical 
dissymmetries of the wind tunnel are exhibited. 
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A resumption of the calculations, but without introducing the hypoth- 
esis of isotropy, produces similar conclusions} the portion of the spec- 
triam (that is, of the spectral tensor) relative to small values of k 
remains constant during the decay, while the rest of the spectrum vanishes 
progressively. It follows that the anisotropy of tvirbiilence that results 
from the geometrical, and mechanical configuration of the flow , and which 
is practically concealed at the beginning of the decay of turbulence, must 
appear in the final phase (Batchelor, ref. 9)* Unfortimately, it is rather 
difficult to check it, since the turbulence is then very weak. 


55. The concept of "dynamic statistical equilibrium" (Heisenberg, ref. 23, 
Batchelor, ref. 8): 

In the problem of the decay of turbulence behind a grill, the study 
of the shapes of spectral c\irves and correlation curves led Heisenberg to 
believe that the notion of similarity could be extended to vortices con- 
taining the major part of the turbulent energy. The overall spectrvim, 
with exception of the region of small wave nimibers or large vortices 
(spectrum in k^), thus would be in a dynamic statistical eq\iilibrlum or 
quasi -equilibri\mi, in which the energy distribution in the spectrum is 
modified d\xring the decay of turbulence, the maximum shifting toward the 
RTTifll 1 wave numbers, but the general shape of the curve remaining the same. 

According to the static theory of local isotropy, the spectral func- 
tion F(k) can be written in the form 


F(k) = Fq^> 



(35-1) 


where — is a length that fixes the transition between the spectral 
ks 

region of predominant inertia forces and that where the viscosity forces 
become comparable to them. <I> is a nondimens ional function. Fq k 

axe dependent on e and v only, and dimensional analysis shows that 

51 _5 1 

Fq = kg = V (35-2) 

It was shown (in section 29) that, when the viscosity is negligible, 

_5 

.F becomes proportional to k which, in this limiting case, fixes 
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the form of the function <5. In the general case, the form of the func- 
tion ii) depends upon that of the function $(k) or T(k) which appears 
in the second term of the fundamental equation (22-8) or (31-1), and which 
represents the transfer of energy through tvirbulence from the large to the 

small size vortices. 3) is in all cases a universal function of — , and 

ks 

the similarity law which it defines is valid in a domain kQ < k < 

where kQ is an . approximate lower limit. This domain is the more extended 


the higher the Reynolds mmiber ^ = — 1-. This number characterizes the 

kQV 

importance of the portion of the spectriim in absolute equilibrium. This 
portion is maintained d\rring the decay of turbulence, so that ^ remains 


constant, as well as the nondimens ional ratio It is found that 

s/^ 

this nimiber intervenes when the region of absolute equilibrium is 
left. The quasi -equilibriimi therefore depends on the parameters e and 

3 2 

V, but, in addition, on a third parameter, the total energy E = ^Uq , of 

which the dissipation e is the derivative with respect to time. By 

Uq^ u a 

(19-1) and (19-2), y is proportional to the Reynolds number R^ = 


Vev 


of turbulence. For a quasi -equllibriimi to be possible, R^ must there- 
fore be constant diiring the development of tiirbulence. This is what the 
experiment proves during the initial phase of decay of turb\ilence, and 
this result supplies an argimient in favor of the quasi -equilibrium hypoth- 
esis. If R^ varies, dynamic similarity ceases, but that does not pre- 
vent the possible existence of a region of local isotropy, variable during 
the development. 


By virtue of the known formula ( 14-11) 


u. 


= 15 


The parameters v 
UqA 


and 


can be replaced by v, A and. u^ 


If 


Sa = 


Uq J. UJ V y /\ W.Q 

is constant in a defined regime of decay, the spectral func- 
tion F(k,t) depends only on a variable parameter, A for example, equi- 
valent to the time, and in addition, proportional to 
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We write 


F(k,t) = F(ki,t)$(n) 


(35-3) 


where r\ = — is a nondimens ional parameter dependent on the fixed wave 
^1 

nmber kj^. The nondimenslonal function is completely defined by the 

value of that is, by the initial conditions at the grid. The length 

of dissipation A can serve as reference length — . The dimensional 

^1 

analysis proves then that F(ki,t) is the product resulting from uq^A, 
or — a quantity proportional to — multiplied by a qmntlty con- 


stant during the decay. Introducing a velocity Vq and a reference 


length — ^ one may put: 

^3 


T] = kA = k 



2 

\2 


F(k,t) = ^1^) $(ti) 


(35-^) 


The form of the function 0 is now defined by means of the fimda- 
mental spectral equation (22-8), which we recall here: 


— + 2vk^F = f 

at 


(35-5) 


We examine first what happens in case of small values of parameter t). 
We shall show that $(ti) is , ^ according to (35-^) ^ infinitely small of the 
first order with respect to t), and that the position of the tangent at 
the origin to the spectral curves F(k,t) is independent of t. For the 
sake of simplification 


F(k,t) = 

\ft 


ni = k\/t 
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By (55-5) 


T)i<l>'l - = 2t^ 


(55-6) 


Assuming that ^ 2 . allows a development of the form 


/ \ n n+1 n+2 


.1 

p 

the first member of (55-6) becomes, except for the factor t , 


(n - 1 )oqT]-^^ + + |(n + l)o 2 + ^ ' 


(55-7) 


£ 

It was seen in section 22, that $(k) was infinitely small as k , 
that is, as . It follows that, either 

n = 1 (X 2 = 0 ( n + 1 ) cx^ + ^ ^ ^ • 

where the first term of (55“T) niust he the term in or else: 


n = 6 


In the hypothesis n = 6, there would res\ilt, for F(k,t) a law k^ 
which is incompatible with the exact law k^ and with experiment. 

For n = 1, there is for F^k^^t^ a law in k which is not compatible 

with experiment when k — >0 hut is verified for small enough values of k, 
except at the limit. In that case 
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— ocqTJj^ + 0^2 + , . . 

and 

F(k,t) = o^k + a^tk^ + . . . (35-8) 


(Xq, 02 being constants independent of t. It follows that 


^ = oq + 3o2tk^ + . . . 


approaches a limit independent of t when k — ^0. All spectral curves 
compatible with the quasi -eqiiilibrium have the same tangent at the origin. 

Therefore, if the part of the spectrum in k^ (represented by 
dotted curves in figure 12) is neglected, the spectral ciorves are approxi- 
mately homothetic with respect to the origin, the maximiam decreasing and 
shifting toward the small values of k when the time increases, that is, 
with progressing decay of turbulence. 

The simdlarlty is obtained again on the correlation curves f and 
g. In fact, (li<--4) 


p» if'i 

R(r) . 2 2iL£!^(k,t)clk = I (55.3) 

Jo rk t Jo _r_ 

" 


The functions Rt, and hence — f, g, depend therefore on the 


^0 

variable and not separately on r and t. In other words, the 

\/^ 

shape of the correlation curves is maintained, provided that or ^ 

Vt ^ 

is taken as variable. 
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Incidentally, ^ can be taken as vaxiable also, because 


L = f dr = 


2Uq2 Jo 


R dr = 


_ rt 


2u, 


0 


'iM a. - ^ rM:ii) 111 

a N/t 


(55-10) 


Since in the region where R^ = C , 
to t, L is proportional to \pt as A, 

But the similarity with L is more quickly destroyed during the decay 
of turbulence than the similarity with A, because in the expression of L, 
the values of corresponding to the small values of are of greater 

weight than the others . The law Uq^A^ = c^® prevails therefore much 

O O "tjG 

longer than the law Uq L = c . The explanation is that similarity is 

approached with respect to L, while it still seems to be verified for A. 
It can be seen that the similarity discrepancies of the correlation cinrves 
appear, first, for the great values of r, and that they afterward reach 
the central zone. The reason for this zone to start changing is that the 

maximum energy of the spectral curve reaches the spectral region in k^. 

The law of decrease of energy can no longer be of the form UQ^t = 

The turbulence reaches the final period of decay, in which the energy u^^t 

_5 

decreases as t 

In this whole analysis it is not necessary to assume that ^ is 
high. The main point is that the initial conditions are such that .the 

region, in which the spectral f\mction F varies as k^, be relatively 
small. If R^ is very high, there exists a region of the spectral curve 

_5 

in which F varies as k 5. This region does not exist when R^ is 
mediimi. 

If it is desired to go farther and define the fonn of the function F, 
the form of the energy transfer f(k) (Heisenberg, reference 2h) must be 


\Iq Is inversely proportional 

and — is constant. 

A 
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chosen. If f(k) has the form given to it by Heisenberg (section 25 ), 
for instance 


T(k) = J(k')dk' 

Jo 



k’^F(k' )dk' 


(55-11) 


the reduced variables ( 35 -^) 3-^® eniployed, and it is proved that ( 35 - 5 ) 
leads to the reduced equation 


iI>(ti' )dTi' 


^n^'(n) 


= 2 


A + C 



Q 


(55-12) 


where A = v— is the inverse of a Reynolds number. 

’^0 

Equation (35-12) is easily reduced to a differential equation which, 
contrary to what occurs in the case of the steadiness hypothesis ( 25 - 6 ), 
seems unsol vable. Heisenberg (reference 2^) has pointed out approximate 
solutions, valid for great or small values of A. (The reader is referred 
to Heisenberg’s report and also to Batchelor, reference 8 ). 


56 . Synthesis of the results relating to the structure of the spectrum 
of turbulence: 

Obviously only isotropic turbulence is involved. The results col- 
lected regarding Heisenberg’s spectral fLinction F(k) which, incidentally, 
cannot be measured directly, but is deduced, by simple transformations, 
from the directly measurable spectral functions (spectral functions of 
Taylor, insofar as they define a spectral distribution in space, and not 
in time ) or from the correlation functions . 

Whatever the structure of turbulence may be, the spectral laws tend, 
when the Reynolds number becomes high enough, toward an absolute limiting 

law, where F(k) is proportional to k 5 (absolute equilibriimi) . The 

2 

corresponding correlation functions are of the form f(r) = 1 - C ^r^. 
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In the problem of the decay of turbulence behind a grid, an initial, 
an intermediate and a final period are distinguished. For small values of 

k, F(k) is always = Ck^, C being a constant independent of time 
(boiind to the invariant of Lo'itsiansky) . 

For the average values of k, and in the initial period, a state of 
"dynamic equilibrium" exists (figure I5) • 

The corresponding portions of the spectraJ. curves are homothetlc with 
respect to the origin and, when the time of decay increases, their maximum 
approaches the origin. Their shape is known up to the origin, but the 
straight part that corresponds to the small values of k and which is 
fixed during the decay of turbulence, has no physical reality and must be 

li 

replaced by the part of the spectrum in k . 


If the Reynolds number is high eno\igh, the portion of the spectrum 
corresponding to great values of k is defined by the laws of absolute 

. . f -I) 

equilibrium \spectrum in k ). When k — ><» the spectrum terminates in 

_5 

a curve k"'^ (experimentally doubtfiil). The region k 5 disappears 
if the Reynolds number is not very high. 


If the Reynolds number is small enough, the region in dynamical 

/ \ 2k^ V 

equilibrium is nimierically known: F(k) is proportional to ke 

p 

provided that k vt is not too great. 


In this initial period the energy of the fluctuations varies as t“^. 


In the final period the zone of dynamic similarity disappears, and 

the meiximum of the spectral cirrve intrudes on the region k^ (fig- 
ure l 4 ) . Starting from an instant t^^ marking the beginning of this 
phase we can write 


F(k) = Ck^e' 2 vk 2 (t-ti) 


aind the energy of the fluctuations decreases as 
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Index of Principal Notations, Fundamental 
Formulas, and Dimensional Equations 

Certain letters have been employed successively, in various paragraphs, 
to represent different quantities without having to be afraid of confusion. 
Others, on the contrary, always represent the same qiiantity; at least from 
chapter II on. We shall enumerate them and recall the principal formulas 
into which they enter and to which one has frequently to refer: 


Xi,X 2 ,X 5 

represent the coordinates of a point M 


designate the differences of the coordinates of two 
points M, M', of which the distance is r. (in 
the first chapter, through a few lines, r repre- 
sents, by exception, a correlation coefficient.) 

f(r),g(r) 

designate Karman's double correlation functions. 

(in chapter I, f is used to represent certain 
statistical frequencies and probability densities . ) 

These functions are derived from the correlation ten- 

in space, of which R = ^aa 

is the scalar invariant. (In the first chapter, R 
denotes the density of probability of position and 
velocity of a statistical point.) R shoiild not be 
confused with R which denotes a Reynolds number 
particularized by a subscript: R^^, ^, . . . 

a(r),b(r),c(r) 

denote the triple correlation functions of Karman 
(section 9) which are deduced from the triple- 
correlation tensor 

Ai,A 2 ,Aj 

are the vector components of the wave number (spatial 
frequency) in the space A of the wave numbers 

k 

is the length of the vector 

F(k) 

is Heisenberg's spectral function, associated with 
the spectral tensor 

E 

denotes the total energy of the fluctxaation of tur- 
bulence per unit mass 


Ikk 
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e denotes the energy dissipated as heat per unit mass 

V denotes the coefficient of kinematic viscosity 

The letter u always designates a velocity. From chapter II on Uq^,U2^u^ 

represent the components of the velocity of the fluctuations at the 
point M, and, if the turbulence is isotropic, uq Is the mean square 

value of u^ . 

X is the dissipation length, L the correlation length 

h always represents a statistical function with orthog- 

onal increments 

The quantities enumerated above have the following dimensional equations: 




F(k) = l 5 t "2 


E = 


l2t-2 


V = 




R = 


2 -2 
L T 


= 


By way of comparison, the dimensional equations of classical mechanical 
quantities read: 

Velocity = 

Acceleration = LT~^ 

Force = MLT ~^ 

Work = ML^T~^ 

Power = ML^T"^ 


NACA TM 1577 


145 


There follows a list of the principal formulas^ referred to con- 
stantly in chapters IV and V. The sections in which these formulas 
appear for the first time are shown in parentheses. 

= ujx)up(x') 


= u, 


0 


(Section 8) < 


f - g 


+ g6„ 


2 0,^3 a3 


^ 1 — ^oux 

a=l 


^a37 = 


(isotropic turbulence) 


A =, -■ 


(Section 9) 


g"(0) 


8f 


L = J f (r)dr 


g = f + (incompressibility) 

2 8r 


(Section I 5 ) 


® = I 


e = 2v k2F(k)dk 
■ 0 


R(r) = 2 


sin r] 


'0 rk 


^(k)dk 


(Section l4)< 


1 P°° 

F(k) = — I rk sin rk R(r)dr 

« Jo 


= 10\£ 
e 
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(Section 19) — = -e 
St 


— Uq^ = -lOv^ 


dt 




— + 2vk^F = 

St 

(Section 22)^ 

$ is a function of k and of t, infinitely small as k^ 

J r) 00 

$ dk = 0 

0 


The following symbols and abbreviations are used: 

X* is the conjugate complex of X. 

u is the mean value of u. Beginning with chapter II the mean 

values refer to a point, considered in the first chapter as 
conditional averages . 

The bold-faced letters in the first paragraph, such as K, represent 
operators . 

For simplifying the writing, notations such as f(u; x) are used in 
place of what should be written unabbreviatedly f ("u 2 _,U 2 ,u^; Xj_,X 2 ^x^^; 

likev/ise, jC‘(u; x) is the abbreviation for 
x^,x^,x^^dUi du^ du^. 


The formulas are numbered in every section; for example, form- 
ula ( 15 -^) is the fourth formula in section I 5 . 


Translated by J. Vanier 
National Advisory Committee 
for Aeronautics 
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APPENDIX 


Some experimental results . 

In the foregoing analytical study of the various recent theories on 
turbulence, reference was made frequently to experimental results, mostly 
of English origin, but numerical values were given only rarely. The data 
contained in this appendix are intended to remedy this to some extent, by 
way of typical examples of several correlation curves and spectrum curves 
obtained by direct observations in the wind tunnel. 

These curves have been determined by A. Favre and his collaborators 
at the Laboratory of Mechanics of the Atmosphere of the Institute of Fluid 
Mechanics, at Marseille, for the official French Aeronautical Research 
Establishment (O. N. E. R. A.). They are copied from reports by 
A. Favre, 13 except the spectral curve .F(k) (fig. 20) which was calcu- 
lated from other data. 

The turbulent velocity was measured with considerably modified hot- 
wire instruments of the Datwyler type which had undergone important modi- 
fications at the Marseille laboratory. 

The time -correlation curve was obtained by means of a recording 
device with Tolana magnetic tape, remodeled in the laboratory so that it 
could be used for aerodynamic purposes. The electric current from the 
hot wire is recorded at a point on the magnetic tape driven at uniform 
forward speed V. Then the record is read at two points separated by a 
distance D, before which the tape unrolls. The currents obtained corre- 
spond to two turbulence records separated by a time displacement ^ and 

then conjugated, as customary in space -correlation measiirement, for the 
two currents coming simultaneously from two distinct hot wires. The 
method also makes it possible to measure the correlations with time and 
space displacement. So, from the experimental point of view, the problem 
of time correlations is solved. In this respect experiment is ahead of 
theory which up to now provides no clear prediction about the time- 
correlation curves . It is to be noted that the correlations at a point 
with- time displacement are more precise to measure than the space corre- 
lations , because only one hot wire is used, and the time displacement 
can be diminished as desired, with the employed instrumentation. By way 
of contrast, for the space -correlations, it is difficult to get two per- 
fectly identical hot wires, and it is impossible to bring them closer 
together beyond a certain limit. 


^5a. Favre: Statistical time -correlation measurements; VII. 

International Congress of Applied Mechanics, London, 1948. 

New time -space correlation measurements downstream from a turbulence 
grid, 0. N. E. R. A. Report No. 8/522-A, Dec. 3 I, I 949 . 
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All the curves reproduced here correspond to an overall speed of 
12.2 m/sec. For the time -correlation curve (fig. l8) the intensity of 


The grid producing the turbulence in figure l6 vas of 5-25 inch 
(8.5cm) mesh width; this value was chosen in order to facilitate the 
comparison with H. L. Dryden’s measurements. In the other graphs a 
grid of 1 inch was involved. In all the cases the spacing of the hot 
wires from the grid was 40 meshes. 

Figure 15 represents a space-correlation curve. The distance r 
of the hot wires, which were placed on a single horizontal perpendicular 
to the tunnel axis, serves as abscissa, Karmans*s dimensionless func- 
tion g(r) as ordinate. 

Figure l6 represents the same function g(r) for a different grid. 
The crosses correspond to a direct measurement; the round points corre- 
spond to the same measurement, but after recording on the time -correlation 
instrument, and rectification with zero time displacement. 

This produces a certain check on the accuracy of the time -correlation 
measurements. The curve was extended up to near 200 mm. It remained 
reasonably constant and below the axis of the r. 

Figure 17 shows the representative curve of the longitudinal corre- 
lation fimction f(r), computed from g(r) by Karman’s formula, and the 

curve representative of the function — (f + 2g), proportional to the func- 

5 

tion R(r), and obtained by calculation from the curve f and the 
curve g of figure l6. The values of R and f are very small and 
negative starting from r = 120 mm and r = 200 mm, respectively. (Com- 
pare formula ( 18-12).) 

Figure l8 represents a time -correlation curve. It is seen to be 
strongly negative for settings above 5 milliseconds. 

Figure 19 represents the spectral function F(n) of turbulence 

obtained by transformation of the time -correlation curves. It is there- 
in 

fore the Fourier transform in cosines of the time -correlation c\irve. 

^Slore precisely, according to Taylor's formulas, F(n) is defined 


time -correlation function, and, in partic\H 


In this case, infinity corresponds practically to 1,000 periods per second. 


turbulence was 1.92 X 10"^. 


so that, if T is the time displacement. 
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It was limited to the low frequencies of around 16 periods per second. 

The theory , as far as it can be applied to the spectral energy distribu- 
tion curves in time at a point, suggests that the spectral curve passes 
through the origin. This corresponds to the fact that the areas of the 
positive and negative parts of the time correlation curve are equal. It 
is verified approximately on figure I8. This verification is not alto- 
gether rigorous, but it does seem that for values of more than 30 milli- 
seconds in time displacement the curve becomes positive again, although 
remaining very close to the time axis . 

Figure 20 gives, computed from the curves of figure I7, the integral 


equal to the quotient of the spectral function F(k) by the constant 5 uq^. 

Therefore, with s\ii table units as ordinates, it is the representative curve 
of F(k)j k is the inverse of a length, while in figure I9, n is the 
inverse of a time Interval. A priori, Taylor's spectral function F(n), 
which represents the spectral energy distribution in time, and Heisenberg's 
spectral fxmction F(k) which represents the spectral energy distribution 
in space, are different. 
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Figure 1.- Recording of the turbulent fluctuation velocity in a wind tunnel 
of 20 cm X 30 cm. Main -flow velocity: 20 m/sec. Intensity of tur- 

bulence: 5-10 . The recording duration is 0.03 sec. (Photograph 

Heubes.) 
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Figure 3 
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Figure 5 
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Figure 9 



Figure 10 
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Figure 11 
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Figure 12 
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Figure 14 
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Figure 15.- Transverse correlation function g(r) in space. 
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Figure 16.- Transverse correlation function g(r) in space. 
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Figure 17.- Longitudinal correlation function f(r). Total correlation 

function f(r) + 2g(r) 
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Fig;ire 18.- Time correlation function. 
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Figure 19.- Taylor’s spectral function F(n). 



Figure 20.- 


Spectral function F(k). ^Value of the integral 
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kr sin kr — dr j 
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